Quantum field theory from classical statistics 
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An Ising-type classical statistical model is shown to describe quantum fermions. For a suitable 
time-evolution law for the probability distribution of the Ising-spins our model describes a quan- 
tum field theory for Dirac spinors in external electromagnetic fields, corresponding to a mean field 
approximation to quantum electrodynamics. All quantum features for the motion of an arbitrary 
number of electrons and positrons, including the characteristic interference effects for two-fermion 
states, are described by the classical statistical model. For one-particle states in the non-relativistic 
approximation we derive the Schrodinger equation for a particle in a potential from the time evolu- 
tion law for the probability distribution of the Ising-spins. Thus all characteristic quantum features, 
as interference in a double slit experiment, tunneling or discrete energy levels for stationary states, 
are derived from a classical statistical ensemble. Concerning the particle-wave-duality of quantum 
mechanics, the discreteness of particles is traced back to the discreteness of occupation numbers 
or Ising-spins, while the continuity of the wave function reflects the continuity of the probability 
distribution for the Ising-spins. 



INTRODUCTION 



(—j Can the quantum-mechanical interference pattern in a 
■4— > double-slit experiment be described by a suitable time evo- 
O lution of a classical statistical ensemble? Can this be ex- 
(1} tended to the quantum- interference of multi-fermion sys- 
terns? Can the evolution law for the classical probabilities 
Be causal and local in the sense that the probability dis- 
I tribution at time t + e can be computed from the one at 
J>time t, and that the evolution of local properties of the 
distribution is only influenced by the distribution in a lo- 
cal neighborhood? We answer these questions with a clear 
''yes" . While a general discussion how "no-go-theorems" 
^ as Bell's inequalities 0, [|| or the Kochen-Specker theorem 
|§] can be circumvented can be found in ref . 0| , we con- 
1— ' centrate in this paper on the construction of a classical 
statistical model that can describe real physical situations. 
. . , We will investigate an Ising-type model for discrete oc- 
^ cupation numbers or Ising-spins on the sites of a lattice. 
KA Intuitively, if a particle is present on a lattice site x the 
j_j Ising-spin is up, and if no particle is present it is down. 
Since the associated occupation numbers only take the val- 
ues one or zero, one sees already the close analogy to oc- 
cupation numbers of fermionic multi-particle systems. The 
notion of presence or absence of particles can be extended 
to more general ground states, as a half-filled state. Now a 
flip of the Ising-spin at x from down to up corresponds to a 
change from a situation with no particle at x to one with a 
particle at x. More precisely, we will consider four or eight 
"species" of Ising-spins. These species correspond to the 
real degrees of freedom for Majorana or Dirac spinors. 

The state of the system or the classical statistical ensem- 
ble is characterized by the probability distribution for the 
possible configurations of Ising-spins. For the formulation 
of dynamics one has to postulate an evolution law which 
describes how this probability distribution changes in time. 
For the Ising-spins there is no notion of underlying continu- 
ous and deterministic classical dynamics for "trajectories" 
-each Ising-spin can only take two values. A differential 



evolution equation has therefore to be formulated on the 
level of probability distributions. This constitutes the fun- 
damental law defining the dynamics. (It is the analogue of 
the Liouville equation for a statistical ensemble of Newto- 
nian classical particles.) The only basic restrictions on the 
form of this law are that all probabilities have to remain 
positive and the sum of all probabilities must equal one 
for all times t. We propose a specific causal and local evo- 
lution law. With this evolution law our Ising-type model 
describes a quantum field theory for Dirac fermions in an 
arbitrary external electromagnetic field. The Schrodinger 
equation for a quantum particle in a potential follows for 
one-particle states in the non-relativistic approximation. 

With our proposed evolution equation for a classical sta- 
tistical ensemble of Ising-spins, the Schrodinger equation 
and all associated quantum phenomena as interference, 
tunneling, the uncertainty relation and discrete energy lev- 
els for stationary states, are implemented within classical 
statistics. We follow here the concept of "probabilistic re- 
alism" that there is one reality, but the fundamental laws 
are formulated as probabilistic laws. It seems not entirely 
excluded that our evolution law for the probability distri- 
bution could be produced by cellular automata [5]. This 
seems, however, rather difficult, and it is not needed for 
our purpose. 

Time evolution laws for classical statistical ensembles 
that account for the full dynamics of quantum particles 
have been found previously for more restricted settings. 
For a single particle in a potential a suitable evolution 
equation for the probability density in phase space has been 
proposed in ref. @. It leads to the Schrodinger equation 
after "coarse graining" . This "quantum evolution equa- 
tion" modifies Liouville's equation. In ref. @ the different 
behavior of classical and quantum particles is traced back 
to the specific form of the evolution equation. A continuous 
interpolation between the "classical" and "quantum" evo- 
lution law can describe "zwitters" - particles whose prop- 
erties interpolate continuously between classical and quan- 
tum particles. This earlier approach constitutes already a 
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powerful demonstration that there is no difference in princi- 
ple between classical statistics and quantum statistics since 
a continuous interpolation between both is possible. It also 
allows the formulation of consistent theories that are arbi- 
trarily close to quantum mechanics, with a small parameter 
accounting for deviations that can be restricted by exper- 
iment. However, the generalization of this setting to a 
multi-particle situation is cumbersome. It is much simpler 
to base a realistic model directly on the analogy between 
Ising-spins and occupation numbers for fermions, as done 
in the present paper. Furthermore, the setting of ref. @ 
does not account explicitly for the discrete particle prop- 
erties, while in the present paper they can be associated 
directly to the discrete Ising-spins. 

In ref. Q, [H, the map between Ising-type models 
and quantum fermions has been developed. In particular, 
ref. Q has already constructed a classical statistical en- 
semble for massless Majorana spinors in four dimensions. 
The present paper builds on these findings. It extends 
this formalism to include a mass term, to formulate Dirac 
spinors and to add the coupling to external electromagnetic 
fields. These steps are needed for a description of realistic 
situations. The discussion of one-particle states can now 
describe electrons in a potential. For a non-relativistic ap- 
proximation we can derive the Schrodinger equation from 
the proposed evolution law for the classical statistical en- 
semble of Ising-spins. 

Our Ising-type model is formulated for a lattice of dis- 
crete points in space and we use discrete time steps. (The 
lattice formulation of the present paper differs in some as- 
pects from ref. 9].) As long as the number of lattice 
points (in space) remains finite all mathematical opera- 
tions are defined unambiguously and the model is fully 
regularized. The map between classical statistics for the 
Ising-type model and quantum fermions can be done on 
this level where no ambiguities are present. At the very 
end one may perform the continuum limits of vanishing 
lattice distance and vanishing size for the time steps. We 
formulate the model in a way such that the unitary time 
evolution is already realized for discrete time steps and a 
finite number of lattice points in space. This implies cer- 
tain restrictions on the precise lattice formulation for which 
more details will be discussed in a separate paper. 

Our paper is organized as follows: In sect. |TT] we for- 
mulate the classical statistical ensemble for Ising-spins and 
we introduce the important concept of the "classical wave 
function" . Up to signs this is the square root of the prob- 
ability distribution. Time evolution laws preserving the 
positivity and normalization of the probability distribution 
can be most easily formulated in terms of the classical wave 
function. Since the normalization condition states that the 
classical wave function is a vector with unit length, we can 
formulate the normalization preserving evolution equations 
as rotations of a vector. We also employ the simple map be- 
tween the classical wave function and a Grassmann wave 
function. It has been shown in refs. Q, @ that a lin- 
ear time evolution of a Grassmann wave function can be 
mapped to a Grassmann functional integral and vice versa. 
This permits us to formulate the evolution law for the clas- 



sical wave function and probability distribution directly in 
terms of the action of a Grassmann functional integral. 
This possibility is of great help for the implementation of 
the symmetries preserved by the evolution as, for example, 
the Lorentz symmetry for relativistic particles. 

In sect. IIIII we present the lattice action for massless 
Majorana fermions in four dimensions. Since the classical 
wave function is a real object, it should be an element of a 
real Grassmann algebra. The continuum limit of the lattice 
action realizes Lorentz-symmetry. In sect. IIVI we extract 
the time evolution of the Grassmann wave function from 
the real Grassmann functional integral defined in sect. IIIII 
It obtains by "integrating" out the Grassmann variables at 
past (or future) times. Sect. IVl shows that the time evolu- 
tion is unitary. This holds already for discrete time steps 
and can therefore be extended in a straightforward way 
to the continuum limit of infinitesimal time steps. Sect. 
fVl establishes the concrete evolution law for the classical 
wave function and therefore the probability distribution 
of the classical ensemble of Ising-spins. In sect. I VII we 
briefly recapitulate how expectation values of observables 
in the classical statistical ensemble can be computed from 
the Grassmann wave function. 

In sect. IVIII we turn to the notion of statistical states 
with a fixed number of fermions. We concentrate on the 
one particle state and show that it describes the prop- 
agation of a massless Majorana or Weyl fermion. This 
can be generalized to multi-fcrmion states. The equiva- 
lence between Majorana and Weyl spinors in four dimen- 
sions gives a first glance on the emergence of a complex 
structure which is discussed in sect. I Villi The presence 
of a complex structure is characteristic for quantum me- 
chanics where the "physics of phases" plays a crucial role. 
While single Majorana spinors are described by a real four- 
component quantum wave function, Weyl spinors are for- 
mulated in terms of an equivalent complex two-component 
quantum wave function. For Dirac spinors a different com- 
plex structure will be used and we therefore resume in sect. 
I Villi briefly the general features of complex structures in 
our setting of a real Grassmann functional integral and a 
real classical wave function. In sect. IIXI we finally add a 
mass term. We present some explicit instructive examples 
for the classical probability distributions that represent the 
propagation of massive Majorana spinors. 

In sect. |X] we turn to Dirac spinors. They are sim- 
ply represented as two Majorana spinors. The electro- 
magnetic gauge transformations rotate the two Majorana 
spinors into each other. The coupling to an external elec- 
tromagnetic field involves therefore both types of Majorana 
spinors. We discuss the complex structure, the usual Weyl 
representation of Dirac spinors and the symmetries of the 
model. From the Grassmann functional integral we ex- 
tract again the evolution equation for the Grassmann wave 
function. The associated evolution equation for the prob- 
ability distribution of the classical statistical ensemble for 
Ising-spins describes the dynamics of an arbitrary number 
of electrons and positrons in an external electromagnetic 
field. Sect. IXII discusses the one-particle states. The clas- 
sical eight-component real wave function obeys the Dirac 
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equation in an external electromagnetic field. The complex 
representation and the non-relativistic limit are straightfor- 
ward. This results in the Schrodinger equation for a par- 
ticle in a potential. All characteristic features of quantum 
mechanics, as interference, tunneling or discrete energy lev- 
els can therefore be derived from our evolution equation for 
a classical statistical ensemble of Ising-spins. This is one 
more clear concrete example that quantum physics can be 
described by classical statistics! Even though this is not 
the main purpose of the present paper, we comment how 
the discrete particle properties in quantum mechanics are 
related to the discrete Ising-spins, and how the continuity 
of the wave function simply reflects the continuity of the 
classical probability distribution. Particle- wave duality ap- 
pears as a simple consequence of our classical statistical 
setting. Our conclusions are presented in sect. IXII1 

II. PROBABILITY DISTRIBUTION AND WAVE 
FUNCTION 

A classical statistical ensemble is specified by its states 
r and a probability distribution {p T }, which associates to 
every r a positive probability p T > 0. The distribution is 
normalized E T p T = 1. Classical observables take in every 
state r a fixed value A T , and the mean value obeys (A) = 
H T A T p T . 

1. Generalized Ising model 

We will consider an Ising-model type system for discrete 
Ising-spins s 7 (x) = ±1, with x points on a suitable three- 
dimensional lattice and 7 denoting different "species" of 
Ising-spins, 7 = 1 . . ,N S , For N s = 4 our system will be 
equivalent to a quantum field theory for Majorana or Weyl 
spinors, while for N s — 8 we will describe Dirac spinors. 
The states r are sequences or configurations of Ising-spins, 
r = {s-y(x)}. Instead of Ising-spins, we will actually use 
occupation numbers or bits n 7 (x) = (s^(x) + l)/2, such 
that the states r describe bit sequences of numbers n 7 (x) = 
0, 1 , t — |n 7 (a?)}. We will consider a setting with L 3 /8 
lattice points (L even) such that the number of states is 
2 (n s l 3 /8)_ The cont i nuum ii m it L — > 00 is taken at the 
end. 

The sequences or configurations of occupation numbers 
{ro 7 (cc)} show already a strong analogy to the basis states 
of a quantum theory for an arbitrary number of fermions, 
formulated in the occupation number basis for positions. 
We will exploit this analogy in order to formulate a funda- 
mental "evolution law" for the time evolution of the prob- 
ability distribution p T (t), such that our system describes a 
relativistic quantum field theory for fermions. 

We emphasize that the form of the evolution law is not 
known a priori - we do not introduce any underlying de- 
terministic theory of "classical trajectories" or similar con- 
cepts. We note that no continuous time evolution for indi- 
vidual Ising-spins s 7 (x) can be formulated, since s 7 (x) only 
admits discrete values ±1. On the other hand, the time 
evolution of the probability density may well be continu- 
ous. We postulate that the basic dynamics describes the 



time evolution of probabilities. The corresponding dynam- 
ical law must be specified in order to define the model. The 
only constraint to be imposed a priori is that is respects 
for all t the positivity and normalization of the classical 
probability distribution 

Pr(*)>0, 5> T (t) = I. (I) 



2. Classical wave function 

A useful concept for the construction of a consistent time 
evolution law is the classical wave function Q , 

q T (t) = S T (t)VpAt) , Pr(t) = q 2 T (t) , S T {t) = ±1. (2) 

This is a real function which is given by the square roots 
of the probabilities up to signs s T (t). Consistent evolution 
laws correspond to rotations of the vector q T (t), 

q T (t) = Y / Rr P (t,t')q p (t'), 

P 

J2 R r P ( t ,t')R*p(t,t') = 6 T(7 . (3) 
p 

The positivity of p T = q% is automatic, and the normal- 
ization of the distribution remains preserved provided it 
was normalized for some initial time ij n , since the length 
of a vector or S T (?^ = S T p T is preserved by rotations. We 
will specify a linear evolution for the wave function. In 
other words, we will consider rotation matrices R Tp that 
are independent of the wave function. 

The classical wave function specifies the probability dis- 
tribution uniquely. The specification of an evolution law 
for the wave function therefore defines the dynamics of the 
classical statistical system completely. In the other direc- 
tion, different choices of the sign functions s T (t) correspond 
to a choice of gauge. For all expectation values and corre- 
lations that can be computed from the probability distri- 
bution {p T (t)} the gauge choice does not matter. Natural 
gauge choices preserve the continuity and differentiability 
properties of the wave function by avoiding arbitrary dis- 
crete jumps which would be induced by arbitrary jumps 
in the sign functions @. In contrast to the discussion of 
classical mechanics in a Hilbert space by Koopman and 
von Neumann [10], the classical wave function is real, such 
that at this step no phases appear as new degrees of free- 
dom beyond the probability distribution. 

Quantum wave functions are usually defined in a com- 
plex Hilbert space. Indeed, many characteristic quantum 
features are closely associated to the "physics of phases". 
In our setting we will define a complex quantum wave func- 
tion by introducing a complex structure in the real space 
spanned by {q T }- The 2^ N " L Z 8 ' real components of the 
vector {q T } can then be associated to a complex vector 
with dimension 2^ NsL Z 8 ) -1 . One complex structure can 
be closely related to the equivalence between Majorana 
and Weyl spinors in four dimensions [Til] , [f2[ [§] , but more 
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general complex structures are possible. There is no con- 
ceptual difference between the classical and quantum wave 
function in our setting. 

3. Grassmann wave function 

We will specify the time evolution of the wave function 
{q T (t)} within a formalism based on a real Grassmann al- 
gebra. This will make the close connection to fermions 
most transparent Q- The Grassmann formulation relies 
on the isomorphism between states r and the basis ele- 
ments g T of a Grassmann algebra that can be constructed 
from the Grassmann variables ip^fa. Each basis element 
g T is a product of Grassmann variables 



^71 0l)V> 7 2 fa) 



(4) 



which is ordered in some convenient way. To be specific, we 
define some linear ordering of the lattice points and place 
variables with "smaller x" to the left, and for each x place 
smaller 7 to the left. If a Grassmann element g T contains 
a given variable ipj(x) we put the number n 7 (x) in the 
sequence r to 0, while we take n 7 (x) = 1 if the variable 
ip-y(x) does not appear in the product ((U). This specifies 
the map between the 2^ NsL / 8 ) independent basis elements 
g T of the Grassmann algebra and the states r. 

An arbitrary element g of the Grassmann element can 
be expanded in terms of the basis elements 
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X<7t3t- 



(5) 



A time dependent wave function {q T (t)} can therefore be 
associated to a time dependent element of the Grassmann 
algebra g(i), provided the coefficients q T (i) are real and 
normalized according to J2 t 1t(^) = 1- We will formulate 
the fundamental evolution law as an evolution law for the 
"Grassmann wave function" git). For this purpose we will 
formulate in the next section a Grassmann functional in- 
tegral for a quantum field theory of massless Majorana or 
Weyl fermions. It will involve N s L 3 /8 Grassmann vari- 
ables i/)^(t,x) for every discrete time point t. In section 
I VI we will present a prescription how the Grassmann wave 
function g(t) and the fundamental evolution law can be ex- 
tracted from this functional integral. This will establish a 
map between a quantum field theory for fermions and Ising 
type classical statistical ensembles with dynamics specified 
by an appropriate evolution law. 



III. MASSLESS MAJORANA FERMIONS 
1. Action 

In this section we formulate the quantum field theory 
of a free Majorana spinor for a discrete space lattice on a 
torus. It will be defined by a Grassmann functional integral 
based on a real Grassmann algebra, and we take N s = 
4. Due to the finite number of (£/2) 3 space points we 
have at any given t a finite number L 3 /2 of Grassmann 
variables ip^(t,x), with four "species" 7 = 1...4. (No 



complex conjugate Grassmann variables are defined at this 
stage.) For the associated classical statistical ensemble the 
states correspond to sequences of L 3 /2 occupation numbers 
n 7 (x) that can take values n^{x) = 0, 1. The total number 
of classical states equals 2 L I" 1 and remains finite for finite 
L. 

We also will take time t on a finite discrete chain. The 
functional integral will therefore involve a finite number of 
Grassmann variables. It is fully regularized and all quan- 
tities are well defined. The continuum limit L — > 00 of an 
infinite number of Grassmann variables will only be taken 
at the end. 

We start with the action for the regularized quantum 
field theory 



s= £ L(t), 



t=Un 

with Lagrangian 

L(t) = ip^(t, x)B 1 (t + e,x). 

X 

Here we define 

B 7 (* + e, x) = Y-ys{{v})Mt + e,*k + «*A), 

{v} 



with 



and 



E=n( e 

M 3 Vj=±i 



(6) 



(7) 



(8) 



(9) 



Y~,s(M) = e 1 - XK + w k I)T k - v lV2 v 3 I . (10) 

Thus B~f{t + e,x) involves a linear combination of Grass- 
mann variables at lattice sites which are diagonal neigh- 
bors of x, corresponding to corners of a cube with basis 
length 2 A, with x being at its center. The sums extend 
over j, k, I = 1 ... 3 and each corner corresponds to a par- 
ticular combination of the three signs vj = ±1. We define 
Wk by wi = — v 2 v 3 ,W2 = viv 3 , w 3 — —viv 2 . Eqs. ©, (fTUj) 
and the following imply a summation over repeated indices 
7,(5 = 1... 4. 

We have also introduced the three real symmetric 4x4 
matrices 71- as 



with 



0,1 
1,0 



To 



0,c 
-c,0 



To 



1 k 



0,1 
-1,0 



= «T 2 , 



T k , 
(11) 

(12) 



and 1 stands for the unit 2x2 matrix. The product of 
these matrices yields the real antisymmetric 4x4 matrix 



c,0 
0,c 



TiT 2 T 3 , I 1 



-I. 



(13) 
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The sum in eq. (|6]) extends over discrete time points 
t n , with f t = ej^t = e Z) n ' — ^« = £,n E 7L, 

tin < £n < £/• The time-continuum limit is taken as 
e — > for fixed ti n ,tf. Similarly, we sum in eq. ([7]) over 
points a; of a lattice. For this purpose we consider for ev- 
ery given t a cubic lattice with lattice distance 2A and 
/ = 8 A 3 We take e — A and place the space-time 
lattice points on a hypercubic bcc lattice with distance 2A 
between nearest neighbors, which we call the "fundamen- 
tal lattice". For even t = 2ne the space lattice points are 
even, x k = 2m' k A, with n,m' k € %. This will be called 
the even sublattice. The odd sublattice consists of the odd 
time points t = (2n + l)e for which the space points are 
also odd, x k = (2m' k + 1)A. The action ([6]) involves indeed 
only Grassmann variables for points on the fundamental 
lattice. For even t eq. (JTJ) involves Grassmann variables 
ip 1 (t,x) living on the even sublattice, while the combina- 
tion B 1 (t + e, x) lives on the odd sublattice. For odd t the 
role of the sublattices is exchanged, now with ip on the odd 
and B on the even sublattice. This construction eliminates 
"lattice doublers" - a more detailed discussion of the lattice 
implementation will be given elsewhere. 

The action ([5]) is an element of a real Grassmann alge- 
bra - all coefficients multiplying the Grassmann variables 
ip 7 (t,x) are real. Within the Grassmann algebra the op- 
eration of transposition amounts to a total reordering of 
all Grassmann variables. The action is antisymmetric 
under this operation, 

S T = -S. (14) 

If we define formally the Minkowski action Sm = iS the 
latter is hermitean, Sm — S M since S M — —Sm- 

As an important ingredient for the probabilistic interpre- 
tation and time evolution discussed in the next two sections 
we observe that one can obtain B~ i from ip^ by a rotation 

v 

^2 ^V( z > y)Rys(x, y) = S„ y S(z, x). (15) 
y 

In eq. (|15[) B 1 and ips refer to the same time arguments, 
but we note that x is not a point on the space-lattice on 
which the Grassmann variables ips(y) live. The quantity 
Bj(x) is a linear combination of Grassmann variables living 
on the corners of a cube with center at x. Since the number 
of lattice points y and the number of centers of cubes x 
is the same, we may nevertheless interpret R^s(x,y) as a 
quadratic matrix. We may formally extend the space to 
objects living on a cubic lattice with lattice distance A for 
each t. In this extended space R acts as a rotation matrix 
despite the fact that x and y refer to points on different 
sublattices of the fundamental lattice. The second equation 
(|15[) and the following equations should be interpreted in 
this sense. 

In order to show that R is on orthogonal matrix we write 
it as a product 

R = R1R2R3, (16) 



with 

R k =D+-T k D k . (17) 
Here the shift operators 

Df(x, y) = ^ [S(x, y-A k )± 6(x, y + A k )] (18) 

act as 

]T Df(x, y)i>{y) = \ yj(x + A k ) ± ^{x - A k )\ . (19) 
y 

They obey 

Y,Dk(*>v) D k(v>*) = ( 2 °) 
y 

i[a(x,z-2A fc ) + 6(x,z + 2A k )±25(x,y)], 

and 

£ D+ (x, y)D- (y, z) = £ D~ (x, y)D+ (y, z) = 

y y 

±[5(x,z-2A k )- 6(x,z + 2A k )]. (21) 

With 

(D±f = ±Df , Rl = D+ + T k D kl (22) 

one finds indeed R k R k = 1- Observing that all shift oper- 
ators , Df' mutually commute, one easily verifies that 
B in eq. (fTQ|) indeed obeys 

B = RiRiR^. (23) 

We next discuss the continuum limits in space and time 
for the action (|5|)- dTUl) . The space continuum limit A — > 
is characterized by 

(D+iP)(x) -> , (D^){x) -> Ad k ^(x), (24) 

with d k = d/dx k the continuum derivative. In the contin- 
uum limit one finds 

R k ^l- AT k d k (25) 

and 

B y (x) = ^(x) - A^2(T k ) yS dMx) +0(A 2 ). (26) 

k 

We further employ 

dMt) = \[i>(t + e)-m], (27) 
such that the Grassmann property tp^ (t, x) = results in 

^~ 1 (t,x)d t ip- y (t,x) = - ip-y(t, x)i> 1 (t + e, x). (28) 
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For A = e we obtain the continuum relation 

i>(t, x)B(t + e,x) = e(ipd t i> - ^ ^T k d k ^). (29) 

fe 

The factor of e combines with E t into J. — eSj. Finally, the 
factor 8A 3 in the relation between f and Y"7 is absorbed 
in the continuum limit by a rescaling of the Grassmann 
variables by a factor (2A)(~ 3 / 2 ). 

2. Lorentz symmetry 

We next show that the action is Lorentz-symmetric 
in the continuum limit e = A, A — > 0. In the continuum 
limit we write the action as 



and obey 



S = / {ipydtipy-tp^iTkj^sdkips}- 



(30) 



It involves now four Grassmann functions ip 7 (t,x),j = 
1 . . .4, x = (x%, X2, £3). The integral extends over three di- 
mensional space and time, with dt = d/dt and d k = d/dx k . 

The Lorentz invariance of the action (|3U|) is most easily 
established by employing the real matrices 



7 



(-" : T i) ^--vn. 



such that 



5 = - / V>7^ , tp = tp 1 i 



(31) 



(32) 



where /i = (0,fc) and <9o = <9 t . We define t/> 7 = i>s(7°)5-y 
(We consider ip here as a vector with components ipj and 
suppress the vector indices. The Pauli matrices are denoted 
by Tfc.) The real 4x4 Dirac matrices 7 M obey the Clifford 
algebra 



(33) 



with signature of the metric given by rj^ v = 
diag(— 1, 1, 1, 1). This can be easily verified by using 
the relations 

{T k ,T l }=2S kl , { 7 °,T fc } = 0. (34) 
Furthermore, one finds 

(7°) T = -7°, (7T=7 fc , (35) 
and the relations 

[T k ,T t ] = 2e klm IT m , [T k ,I] = , P = -1, 
7 o 7 i 7 2 7 3 = J ; ^ J } = ) {<y * J } = Q (36) 

Infinitesimal Lorentz-transformations can be written as 

5ip-y = ~^ e M^ ( s ' il/ ) 7 5 , = -e^, (37) 

where we have omitted to indicate the associated transfor- 
mations of the coordinates. The Lorentz generators 
obtain from the Dirac matrices as 



E /- = _l[ 7 /^] 



^Ofe 



1 

, 



Ti \^kl kirn jrri 

k 1 = — -6 ii„ 



(39) 



We recognize in eq. (1321) the standard Lorentz invariant 
action for free Majorana fermions in a Majorana represen- 
tation of the Clifford algebra with real 7 M -matrices. 

3. Functional integral 

The functional integral is defined by the partition func- 
tion 



Ttygf [tp(tf)]e- s g in [ip{t in )], 



with the functional measure 



Vif> = JJ / (dtp 4 {t,x) .. .dipi(t,x)) . 



(38) 



(40) 



(41) 



The boundary terms gi n and g f only depend on the Grass- 
mann variables ipiUn) and ip(tf), respectively. As we will 
see below, the boundary terms gj and gi n are related to 
each other, such that the functional integral (14TJ1) is fully 
specified by the choice of gi n . 

The basic definition of Z is formulated on the discrete 
space-time lattice with finite volume. Thus Z is a well 
defined real number. We will show later that for a suitable 
normalization of gi n one obtains Z — 1, independently of 
L. The continuum limit L — > 00 for the partition function 
is therefore trivial. 

We will use the Grassmann function integral (I4U|) in or- 
der to specify the fundamental evolution law for the proba- 
bility distribution {p T (t)} of the classical statistical ensem- 
ble for Ising-spins. The positivity and normalization of the 
probabilities holds for an arbitrary choice of gi n , provided 
g~f is related to gi n appropriately. For every given gi n the 
probability distribution {p T {t)} is uniquely computable for 
all t, such that the functional integral (14T)1) indeed specifies 
the time evolution of the probability distribution. 



IV. GRASSMANN WAVE FUNCTION FROM 
FUNCTIONAL INTEGRAL 

In this section we compute for the functional integral 
(I40p a Grassmann wave function g(t), which is an clement 
of the Grassmann algebra constructed from the Grassmann 
variables vp 7 (t,x) at given t. The central idea is to "inte- 
grate out" the Grassmann variables at times t' ^ t [i"3| . 
The expansion coefficients of g{t) will specify the classical 
wave function {g T (£)}. 

1. Integrating out past and future 

The functional integral (l40l) involves variables for arbi- 
trary time points t n . In order to construct a wave function 
g(t) which only refers to a particular time t we have to "in- 
tegrate out" the information referring to other time points 
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f^t Jim. This can be done by decomposing the action 

S = S<^ ~\- »S'> ; 

S < = ^L(t>) , 5> = 5^L(0. (42) 
t'<t t'>t 

The wave function g(t) obtains now by integrating out all 
Grassmann variables for t' < t 

g(t) = J Vi>{t' < t)e- s <g m . (43) 

We observe that g(t) depends only on the Grassmann vari- 
ables ip{t). More precisely, it is an element of the Grass- 
mann algebra that can be constructed from the Grassmann 
variables ^> 7 (i, x). 

The conjugate wave function is defined as 

g(t) = J VW > t)g f e- s > . (44) 

Again, this is an element of the Grassmann algebra con- 
structed from il>(t). In terms of g and g the partition func- 
tion reads 

Z = J Vil>(t)g(t)g{t), (45) 

where the Grassmann integration J Dif>(t) extends now 
only over the Grassmann variables ij)y(t, x) for a given time 
t. 

2. Classical probabilities and wave function 

We next expand g in terms of the basis elements g T of 
the Grassmann algebra generated by the variables ipy(t), 

g(t) = Y l ( ir(t)gr[m]- (46) 

T 

We associate the real coefficients q T (t) with the classical 
wave function, such that the classical probabilities obtain 
as p T (t) = Qr(t)- This requires for every t the normalization 
^ T <Z?(i) ~ 1- We will show in the next section that this 
normalization condition is indeed obeyed, provided it holds 
for the initial wave function g(ti n ) = gi n [i/) 7 (^„, x)l . 

The conjugate basis elements of the Grassmann algebra 
g T are defined [7( by the relation 

g T gr = niI^W ( 47 ) 

x 7 

(no sum over r) and the requirement that no variable ip-y(x) 
appears both in g T and g T . They obey 

J Th/>(t)g T [m] 9p [m] = Sr P . (48) 

Expanding 

g(t) = J2^(t)9r[m] (49) 



yields 

Z = J2Qr(t)q r (t). (50) 

T 

We will see in the next section that for a suitable choice of 
gf the conjugate wave function obeys for all t the relation 
q T {t) — q T (t). Together with the normalization condition 
J2 T Qr = 1 this guarantees Z = 1. 

In consequence, we can express the classical probabilities 
p T (t) directly in terms of the Grassmann functional integral 

Pr(t) = J V^g f V T {t)e- 3 g in , (51) 

with V T (t) a projection operator 

V T {t)g p [i>{t)] = g T [^{t)]8 Tp , 

J V^(t)g a [i>(t)]'Pr(t)g P [m] = SraSrp- (52) 

The formal expression of V T in terms of the Grassmann 
variables ip(t) and derivatives d/dip(t) can be found in ref. 

Q. 

The interpretation of the Grassmann functional inte- 
gral in terms of classical probabilities is based on the 
map g(t) — > |p T (i)}, which in turn is related to the map 
{qA -> {pA = {q 2 A- The map g(t) O {q T {t)} is invert- 
ible. A given state or classical statistical ensemble may be 
specified by the "initial value" at some time to,g(to), or the 
associated wave function {q T (to)} or classical probability 
distribution {p T (to)}- This is equivalent to the specifica- 
tion of gi n = g(U n ) in the functional integral. 

V. TIME EVOLUTION 

In this section we compute the time evolution of the 
wave function {q r (£)} and the associated probability dis- 
tribution {p T (t)} — {<j^(t)}. This will lead to a type of 
generalized Schrodinger equation for the real wave func- 
tion {<7 r (i)}, as well as an associated evolution equation 
for the Grassmann wave function g(t). We will use the 
properties of this evolution equation in order to establish 
that the norm ^2 T q^.(t) is conserved, such that {q^(t)} 
can indeed be interpreted as a time dependent probability 
distribution. 

Due to the particular form of the action ([6]) , which only 
involves one type of Grassmann variables (and no conju- 
gate variables as in ref. [7j) the discrete formulation of the 
functional integral is crucial. We will see that g(t) jumps 
between two neighboring time points, while it is smooth 
between g(t+2e) and g(i). We will therefore distinguish be- 
tween even and odd time points and use the definition (|46l) 
of the wave function only for even times. (For odd times 
one may employ a different definition, which will guarantee 
the smoothness of the time evolution of {q T (t)} for both 
even and odd time points [9].) In the continuum limit, 
e — > 0, the time evolution of the wave function is described 
by a continuous rotation of the vector |g T (£)}. 
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As a consequence of its definition (l46l) the Grassmann 
wave function obeys the time evolution 

g(t + e) = J V4,(t)e- L(t) 9(t). (53) 

This determines {q r (t + e)} in terms of {q T (t)}. Thus the 
action (jrJj)- (fTU)) specifies the dynamics how the probability 
distribution {p r (i)} evolves in time. The particular dy- 
namics of a given model is determined by the form of B~ f 
in eq. (flOj). 

1. Unitary time evolution 

We want to establish the relation <j3j) which guaran- 
tees the preservation of the sum of probabilities p T = 
J2 T It — 1; wrtn t' t, t —> t + 2e. For this purpose we 
employ the property ([15]) . The relation between q T (t + 2e) 
and q T (t) is computed from cq. (|53[) . 

Inserting the specific form ([7]) for L(t) the evolution 
equation (l53l) reads 

ff(t + e) = /" PV>(*)exp{-^^^M£ T (*+ e ^)}s(t). 

J x 7 

(54) 

We may write <?(£) in a product form 

a(t) =Y[Y[[a-/(t,x)+b 7 (t,x)^(t,x)], (55) 

with some fixed ordering convention of the factors assumed, 
e.g. smaller 7 to the left for given x, and some linear 
ordering of the lattice points, with "lower" points to the 
left. In the product form eq. (l54l) yields 

9(t + e)= fvmilH (56) 

^ x 7 

I [l — ipj(t, x)B 1 {t + e, a;)] [a 7 (i, a;) + & 7 (t, x)tp^(t, x)] j 
= IT IT {^t(*> - t ) + Vj a y(t^ x)B 1 (t + e, ac)|. 

x 7 

Here we use the fact that each individual Grassmann inte- 
gration f dip-y(t,x) can be performed easily, 

J dip(l - ip(f)(a + bip) = b- aip, (57) 

and ?7 7 = ±1 results from the anticommuting properties of 
the Grassmann variables ip, with 

Vi = ? /3 = 1 , V2 = = -1- (58) 
As a result, we can write 

9{t + e) = <?r (*)C.9t [5 7 (* + e, &)] , (59) 

T 

where C<? r obtains from g T by the following replacements: 

(i) for every factor ip 1 {x) in g T one has a factor 1 in Cg T ; 

(ii) for every pair (a;, 7) for which no ipj(x) is present in g T 



one inserts a factor T)yip~(x) in C(/ T . The ordering of the 
factors rj-fip^x) is the same as the ordering assumed in the 
product (|55|) . This implies that we can indeed write the 
action of C on the product (|55|) as 

esc*) = n n (m*> s) + »w*, a# 7 & x )) ■ ( 6 °) 

a: 7 

The conjugation operator C maps every basis element g T 
into its conjugate element g T up to a sign <7 r = ±1, 

Cg T = <Jt9t- (61) 

Applying C twice on the product (|55[) multiplies each factor 
by rj-y. The factors ?y 7 drop out due to the even number of 
minus signs, such that C is an involution, C 2 = 1, or 

C 2 g T =g T . (62) 

The jump between g and Cg for neighboring time points 
suggests the use of eq. (|46]l for the definition of the wave 
function {q T (t)} only for "even time points", namely those 
that obey t n — £j„ + 2ne, n € N. Repeating the procedure 
leading to eq. (|59|) one obtains 

g{t + 2e) = Y / Qr (t)ffr [A (t + 2e, x)] , (63) 

T 

where ^4 7 is a linear combination of Grassmann variables 
tps(t + 2e,x) which live on the same sublattice of the fun- 
damental lattice as ips(t,x). For an arbitrary linear trans- 
formation 

B 7 (t + e,x) = F jS (x,y;t + e)ijs(t + e,y) (64) 
with unit Jacobian, det-F — 1, we can write 

g(t + 2e) = J V^(t + e) 

exp { — ifjry (t + e, x)B 1 (t + 2e, x)\g(t + e) 
= J VB(t + e) (65) 
exp { - B y (t + e, a;)A 7 (i + 2e, x)}g(i + e), 

with 

,4 7 (t + 2e,x) = {F- l )T 5 {x,y;t + e) 

xF Sri {y,z;t + 2e)%(t + 2e,z). (66) 

Here we sum over repeated indices S, rj and repeated co- 
ordinates y, z and write the Grassmann integration as an 
integration over new variables B 7 (t + e,x). We next use 
the expression (|56p for g(t + e) and perform the integration 
over B 1 (t + e, x), 

g(t + 2e) = 1 [[T] 1 (a 1 (t,x) + b 1 (t,x)A 1 (t + 2e,x)). (67) 

x,1 

This replaces in g(t) every factor ip 7 (t, x) by A 1 {t + 2e, x) 
such that eq. specifies the expression A 7 (t + 2e,x) 

appearing in cq. (1531 . 
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We next employ the important property that F is given 
by the orthogonal matrix R in eq. (|15l) . which is inde- 
pendent of t. Since (F~ 1 ) T — F we end with the simple 
expression 

A 1 {t + 2e,x) =B% r) (x,z)i/> v {t + 2e,z), (68) 

where 

R 2 n (x,z) = R jS (x,y)Rsr,(y,z) (69) 

defines a transformation between Grassmann variables on 
the same sublattice. The matrix R 2 is, in turn, also or- 
thogonal (R 2 ) T R 2 = 1. We conclude that g{t + 2e) can be 
obtained from g(t) by a simple rotation of the Grassmann 
variables. 

Since A 7 (x) is related to ip 7 (x) by a rotation it is 
straightforward to show that g T [A 7 (a;)] is also connected 
to g T [tp 1 (#)] by a rotation among the basis elements 

g T [A y (x)] = ^g p [^ 7 (x)]i? pT , 
p 

^^RrpRap = S Ta . (70) 

P 

One infers 

g(t + 2e) = J2^{t)9p[^{t + 2e)]R pr 

= Y / 1r{t + 2e)g T [^ 1 {t + 2e)}, (71) 

T 

with a rotated wave function 

q T {t + 2e) =Y / Rr P q P {t). (72) 
p 

This establishes eq. ©. 

Rotations preserve the length of the vector {q T } such 
that J^tIt^) 1S independent of t. Choosing g in = 

J2r1-r(tin)9r[lp(tin)] with J2r Qrfan) = 1 one infers 

£ T g£(t) = 1 for all t. Therefore {p T (t)} = {q 2 (t)} has 
indeed for all t the properties of a probability distribu- 
tion, namely positivity of all p T and the normalization 
p T = 1. In analogy to quantum mechanics we call a 
time evolution which preserves the norm of {q T (t)} a "uni- 
tary time evolution" . A unitary time evolution is crucial 
for the probabilistic interpretation of the functional inte- 
gral pOl) . 

2. Evolution of conjugate wave function 

We next want to show the relation (for t even) 

g(t-2e)=J2Qr(t)Rr P g P ['>Py(t-2e,x)]. (73) 

r, P 

The definition (T4~4l of the conjugate wave function q then 
implies 

q T (t-2e)=J2q p (t)R pr , (74) 
p 



such that 

q T (t)=Y,Rr P q P ( t -^)- (75) 

P 

Comparing with eq. (1721) one infers that q T (t) and q T (t) 
obey the same evolution equation. (By an analogous argu- 
ment eq. (fT5"j) also holds for t odd.) If q and q are equal for 
some particular time to, they will remain equal for all t. 

If {q T (to)} equals {q T (to)} for some time to we can use 
q T (t) = q T (t) for all £ and infer from eq. (|50l) 

Z = (76) 

r 

As it should be, Z remains invariant under rotations (|T2[) 
of the vector {q T } and is therefore independent of t. 

In order to show eq. (1731 we employ the definition (|44[) 
which implies 

~9(t-e) = Jv^(t)g(t)e- L{t - e \ (77) 

or 

g(t - e) = / Vyj(t)g(t) exp { - £ ^ 7 (t - e, i)B 7 (t, i)} 

= / 2? tp{t)g{t) exp { -^^^-^^■K^y)^^,?/)} 

- fvip(t)g{t)exp{-Y,Bi(t-e,x)1iy(t,x)}, (78) 
with 

B 7 (t,x) = ^2^ s {t,y)R s .y(y,x) 

y 

= ^2(R-% s {x,y)Mt,v)- (79) 
Performing the Grassmann integral one finds 

s(* - e ) = [^(* - e > *)] . ( 8 °) 

r 

where the map C acts similarly as C, with ?7 7 replaced by 

V~i = -VjtVi = m = -1>% = % = 1) and C 2 — 1. A 
similar step yields 

§(t - 2e) = ^ 9r(t)5r [A\(t - 2e, x)] (81) 

T 

with 

A^t - 2e, x) = (i? 2 )^(x, j/)^(t - 2e, y). (82) 
One concludes 

g T [A\ (t-2e,x)]=J29 P (* " 2e, *)] #p~t (83) 
p 

and infers eq. (|T3"j) . 
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3. Boundary terms 

The final point we have to settle in order to establish the 
normalization Z — 1 and the expression (|5ip concerns the 
equality of q(to) and q(to) for some arbitrary time t . This 
is achieved by a proper choice of the relation between the 
boundary terms g~f and gi n in eq. (|40[) . For this purpose we 
may imagine that we (formally) solve the evolution equa- 
tion ((72)) in order to compute \q T (tf) } in terms of { q T (ij„ } , 

9in = J2t qr(tin)gT [ip(tin)] ■ (We assume that t in and tf 
are even.) Using 

9(tf) = ^qr(t f )gr[iKtf)], (84) 

T 

9(tf) = X)9r(*/)flr[^(*/)] =§/, (85) 

T 

it is sufficient to choose g~f such that q T {tf) — q T (tf). 
Equivalently, we may specify the wave function { q T (t ) } = 

{<?t(^o)} at some arbitrary even time to and compute the 
corresponding gi n and t// by a solution of the evolution 
equation, using the fact that the rotation ([72]) can be in- 
verted in order to compute q[t — e) form q(t). 

If we would not adopt the choice q T (tf) — q T (tf) the 
functional integral (|4Tfl) would amount to a transition am- 
plitude, which is another useful notion in quantum field 
theory. We are interested, however, in a classical proba- 
bilistic setting and therefore focus on the choice of g~f that 
guarantees Z = 1 for all t. In practice, neither gi n nor 
g~f need to be computed explicitly since we only use the 
Grassmann functional integral for extracting the evolution 
equation for the Grassmann wave function g(r) and the 
classical wave function {q T (t)}. 

4. Continuous evolution equation 

Finally, we cast the evolution law (1721) into the form of a 
differential time evolution equation by taking the limit e — > 
0. This results in a generalized Schrddinger type equation 
for the real wave function |g T (i)}, 

d t q T {t)=Y J K Tp q P (t). (86) 
p 

Since the evolution describes a rotation, the matrix K is 
antisymmetric 

K pT = -K Tp . (87) 

We identify this evolution equation with the generalized 
Schrddinger equation for a quantum wave function for the 
special case of a real wave function and purely imaginary 
and hermitean Hamiltonian H = iHK. 

The time evolution (156")) translates directly to the prob- 
abilities (no summation over r here) 

d t p T = 2 ^ KrpSrSpy/PrPp- (88) 

p 

Once the signs s T (to) are fixed by some appropriate con- 
vention at a given time to, the signs s T (t) are computable 



in terms of the probabilities p T (to)- This follows since for 
all t the wave function q T (t) is uniquely fixed by q T (to) or 
p T (to), andp r (i) is uniquely determined by q T (t). In princi- 
ple, it is therefore possible to formulate the time evolution 
law for the probabilities uniquely in terms of the proba- 
bilities p T (to). However, an expression of dtp T (t) in terms 
of p T {t) needs to keep track of the sign functions s T (t). 
In principle, this can be done by an updating procedure. 
The sign s T can flip only at times t where p T (t) = 0. If it is 
flipped or not at these times is decided by the requirements 
p T > 0, S r p T = 1 for the following times [9j. 

Instead of such an updating procedure it is much more 
convenient to use the wave function and the linear evolu- 
tion law ([86]) for the description of the classical statistical 
ensembles associated to the action S. The basic reason is 
the condition of unit norm of the probability distribution. 
This can be quite cumbersome for a general evolution equa- 
tion for {p T }, but it is extremely simple on the level of {q T } 
where only the length of a real vector has to be preserved. 

5. Grassmann evolution equation 

The matrix K can be extracted from the Grassmann 
evolution equation 

d f g = Kg, (89) 

according to 

dtg(t) = ±- e [g(t + 2e)-g(t)]=ICg(t) = Y / <lAt)ICg T 

T 

= 9r (t)g P ( x )] k pt = ^2d t q T (t)g T [ip-y (x)] . 

(90) 

(In eq. (|9"0)) we use a fixed basis, corresponding to the 
basis elements g T constructed from tjj(t) for g(t), and from 
tp(t + 2e) for g(t + 2e).) For our model of free fermions the 
Grassmann evolution generator K, is given by 

This yields the matrix element in eq. (|86|). 

Kpr = J Vyjg p K,gr . (92) 

In order to proof the relations (|M1) , (liTTj) we infer from 
eqs. (|63|) . (|68|) the relation 

g(t + 2e) - g(t) = £ q T (t)(g T [R 2 ^] - g T &}) (93) 

T 

and use the continuum limit e = A — > 0, cf. eqs. (|16p . 

(USD, 

R 2 ii=(l-2eJ2Tkd k )i>- (94) 
k 
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We then employ the identity (in linear order in e and sum- 
mation over k implied) 



g T [(l - 2eT k d k )iP] 



( 1 + 2e ?*4 
(l + 2e/C)ft#] 



(T k )^ s d k ip s (x) g T [ip] 



(95) 



in order to establish eqs. (159")) , (|9"Tj) . 

Taking finally the space-continuum limit by rescaling 
ip-y(x) and d/dip^(x) such that 



5 



dip 1 {x) 



i>s(y) > = S lS S 3 {x - y), 



(96) 



we arrive at the continuum form of the Grassmann evolu- 
tion equation 



dtg = ICg , JC = 







dtjjj (x) 



(T k ) jS d k Mx)- (97) 



This evolution equation will be the basis for the interpre- 
tation of the time dependent wave function {q T {t)} and 
probability distribution {pr(i)} in terms of propagating 
fcrmionic particles. 



This yields the quantum rule for expectation values 

(A) = (qAq) =J2<lrA Tp q pi (101) 

with A a diagonal operator A rp = A T S Tp . In the Grass- 
mann formulation one uses the associated Grassmann op- 
erator A obeying 



such that 



Ag T = A T g T , 



(A) = J V^gAg. 



(102) 



(103) 



Here g is conjugate to g, i.e. for g — ^2 T q T g T one has 
9 = E r 2t<?t- 

In classical statistics the time evolution of the expecta- 
tion value is induced by the time evolution of the proba- 
bility distribution 



(A(t)) = "£p T (t)A T . 



(104) 



This corresponds to the Schrodinger picture in quantum 
mechanics 



VI. OBSERVABLES 

Classical observables A take a fixed value A T for every 
classical state r. In classical statistics the possible out- 
comes of measurements of A correspond to the spectrum 
of possible values A T . The expectation value of A obeys 



{A) =Y,PrA T . 



(98) 



Our description of the system will be based on these clas- 
sical statistical rules. For example, we may consider the 
observable measuring the occupation number N 7 (x) of the 
bit 7 located at x. The spectrum of possible outcomes of 
measurements consists of values 1 or 0, depending if a given 
state r = [n 7 (x)] has this particular bit occupied or empty. 

For the Grassmann basis element g T associated to r 
one finds N 7 (x) = if g T contains a factor ip^lx), and 
N 7 (x) = 1 otherwise. We can associate to this observable 
a Grassmann operator A/" 7 (x) obeying (no summation here) 

d 

N 1 {x)g T = (N^(x)) T g T , JV 7 (x) = — ^ i(x)- (99) 

Two occupation number operators Af~ (1 (xi) and N l2 {x2) 
commute. 

In general, we may associate to each classical observ- 
able A a diagonal quantum operator A acting on the wave 
function, defined by 



(Aq) T = A T q T . 



(100) 



(A(t)) = (q(t)Aq(t)), 



(105) 



or the corresponding expression in terms of the Grassmann 
algebra 



(A(t)) = / Vi>~9{t)Ag{t). 



(106) 



Using d t q = Kq (|86|) and d t g = ICg (|97| we infer for the 
time evolution of the expectation value the relations 



d t (A) = (q[A ) K}q)= / V^g[A,JC]g 



(107) 



Conserved quantities are represented by Grassmann oper- 
ators that commute with /C, [A, IC] = 0. 



VII. PARTICLE STATES 

Our system admits a conserved particle number, corre- 
sponding to the Grassmann operator J\f, 



d 



y dtpj{y) 



1>Jy) , W,JC}=0. (108) 



The particle number is Lorentz invariant. We can decom- 
pose an arbitrary Grassmann element into eigenstates of 
N 



g = ^ A m g m , Ng m = mg n 



(109) 
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The time evolution does not mix sectors with different par- 
ticle number to, such that the coefficients A m are time in- 
dependent 



d t g = ^2 A m d t g m , d t g m = Kg v 



(110) 



We can restrict our discussion to eigenstates of N ' ■ The 
range of to is [0, B], with B — N s L 3 /8 = L 3 /2 the number 
of independent Grassmann variables. 

1. Vacuum 

Let us consider some static vacuum state go with a fixed 
particle number mo, 

Kg = , Afg = m g , / Vijjg g = 1. (Ill) 



An example for a possible vacuum state is the totally empty 
state go = |0), with 



io>=n^ = nn^(*)=n^ 2 ^)' 



x 7 



AA|0) = 0. 



It obeys 



xty|o> = i , |6) = i. 



(112) 



(113) 



For too 7^ we shift the particle number by an additive 
"renormalization" n = to — to , such that the vacuum cor- 
responds to n = 0, and g = A n g n . An eigenstate of A/" with 
eigenvalue to = mo + n is called a n-particle state if n is 
positive, and a n-hole state for negative n. 

2. One-particle and one-hole states 

We next define creation and and annihilation operators 



ajj,(x), a 7 (x) as 
aUx)g 



d 



dipj(x) 



g , a 7 {x)g = ipry(x)g. 



(114) 



They obey the standard (anti-)commutation relations 
{a\(x), a t (y)} = 8 ie 5(x - y) , Af = at / (x)a J (x), 



[a\{x),AJ] = -a\(x) , [a 7 (x),N] = a 7 (x). (115) 

Acting with the creation operator on the vacuum produces 
one-particle states 

gi(t)= / q~/(t,x)at / (x)g = Sigo, 

J X 

(N - m )gi = g x . (116) 

If needed, we may multiply Q\ with an appropriate nor- 
malization factor such that the wave function obeys 



(117) 



and g\ has a standard normalization. 

Similarly, a one- hole state with n = — 1 obtains by em- 
ploying the annihilation operator 



g-i(t) = / q~ t (t,x)a 7 (x)g 

J X 

= q 1 {t,x){-i°) 1 sas{x)go^Q-igo, 



(J\f -m )g-i 



-g-\- 



(118) 



No one- hole states exist for the vacuum (|112l) . but this 
issue is different if mo =^ 0, as for example for go — 1 where 
to = B. 

If we transform the one-particle wave function q 1 {t,x) 
infinitcsimally according to 



(119) 



the operator Q\ is Lorentz invariant. (We omit here the 
part resulting from the change of coordinates.) Similarly, 
q~y transforms as q 7 and the corresponding infinitesimal 
transformation results in an invariant G—\. If the vac- 
uum go is Lorentz invariant, the Lorentz transformed one- 
particle or one-hole wave functions will obey the same evo- 
lution equations as the original wave functions. 

The time evolution of the one particle wave function g 7 
is given by 



d t gi = ICgi = J (d t q-^-)g = J , 



r d 



5o- (120) 



and similar for the hole. With 
d 



dtjj 7 (x) 



dk 



dlp e {x) 



(?fe) £7 



and 



[K.,ip 7 (x)] = -d k ip e (x){T k ) 



(121) 



(122) 



one obtains Dirac equations for real wave functions (with 
d = d t ) 



Yd^q = , j^q = 0. 



(123) 



We emphasize that these equations follow for arbitrary 
static states go which obey K,go = 0. 

3. Weyl and Majorana spinors 

Let us consider go = |0) where a 7 (a;)|0) = implies that 
no hole states exist. There are then only particle states 
and the propagating degrees of freedom correspond to Ma- 
jorana fermions. In four dimensions Majorana spinors are 
equivalent to Weyl spinors [ll|. Indeed, we may introduce 
a complex structure by defining a two-component complex 
spinor 



ip{x) 



(px{x) 
<p 2 (x) 



ipi = qi+ iq 2 , <^2 = <73 + «<?4- (124) 
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The matrices T k = j°-f k are compatible with this com- 
plex structure. They are translated to the complex Pauli 
matrices 



T, h 

t = 77 -> Tk, 



(125) 



while the operation (p — > iip corresponds to the matrix 
multiplication q — > Iq. (Cf. ref. [9j for more details of the 
map from Majorana to Weyl spinors.) The Dirac equation 
reads in the complex basis 

d t cp = r k d k tp , id t if = -T k V k ip , V k ip = -id k ip, (126) 

and the Lorentz generators are given by 

(127) 



^ 2 ' 2 

In contrast, the multiplications with 7 cannot be repre- 
sented by a multiplication of ip with a complex 2x2 matrix, 
since 



7 9 = ¥>■ 



-T 2 <p . 



(128) 



If we express ip* in terms of the two-component complex 
vector 



X = E<p = -zr 2 ip =^ 
the transformation q — > 7°q corresponds to 



(129) 



(130) 



We may now introduce the complex four component vec- 



tor 



M 



(131) 



for which all matrix multiplications q — > 7^5 can be rep- 
resented by multiplication with complex matrices of the 
Clifford algebra, 



7 



0, -i 
-i, 



> 7 



0, -ir k 
irk, 



(132) 



1 ( n, 



2 V 0, -T k 

We can also define the matrix 



^fcZ _ _ l _ e klm ( T m, 



0, T„ 



7 



• 1 2 3 

-27 7 7 7 



1, 
0, -1 



(133) 



for which ip and \ are eigenvectors with eigenvalues ±1. 
(Often 7 is denoted as 7 5 .) The transformation q — > Iq 
acts on as ipM — > iyipM- Thus the representation (|133() . 
7 = —il, is consistent with the complex structure. Since x 
is not independent of ip the spinor ^ m obeys the Majorana 
constraint fill]. [Tl. 



B-^* M = * M 



B = B~ 












1 








-1 








-1 








1 












(134) 



where Bj^B^ 1 = (7^)*. 

The real Dirac matrices J( M \ m the "Majorana repre- 
sentation" (|31[) and the Dirac matrices l^ w ^ m the "Weyl 
representation" (|132|) are related by a similarity transfor- 
mation, 



1 



<l : : — A^M , 7(V) = ^(M)^ ( 135 ) 



with a unitary matrix A 

/ 1, 0, ,0, 1 

V 0, -t, -t, 0. 



4. Weyl particles 

The propagating degrees of freedom correspond to the 
solutions of the evolution equation. They are most simply 
discussed in terms of the complex equation (|126p . We can 
perform a Fourier transform 



<p(t,x)= / <p(t,p)e*»- 



^L^(t,p) e <*», (137) 



such that the evolution equation becomes diagonal in mo- 
mentum space, 

id t ip(t,p) = -p k T k ip{t,p). (138) 

The general solution of eq. (|126p obeys 

<p(t,p) =exp(ip k Tkt)tp(p), (139) 

with arbitrary complex two-component vectors <p(p). This 
is the standard time evolution for a propagating Weyl par- 
ticle. 

5. Multi-fermion states 

States with arbitrary n describe systems of n fermions. 
This is not surprising in view of our translation of the clas- 
sical statistical ensemble to a Grassmann functional inte- 
gral. For go = |0) there are only n-particle states, and no 
hole states. They can be constructed by applying n cre- 
ation operators alix) on the vacuum. For example, the 
two-fermion state obeys 



92(t) = -7= / Q-ye(t,x,y)a\(x)al(y)g . (140) 



Due to the anticommutation relation 

d d 



{a\(x) , at(y)} = 



d^(x) ' dip e (y) 



(141) 



the two-particle wave function is antisymmetric, as appro- 
priate for fermions 



q 7e {t,x,y) = -q e7 (t,y,x). 



(142) 
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The Grassmann elements 52 describe states with two Weyl 
spinors. A particular class of states can be constructed 
as products of appropriately normalized one particle states 

g(D,g(a) 

q le (x,y) = q^{x)qf\y) - q^(y)q^(x). (143) 

General two particle states are superpositions of such 
states. 

6. Symmetries 

Besides Lorentz symmetry the action (1301) is also invari- 
ant under global SO(2) rotations 

1^1 — cos a ipi — sin a ip2 , ip'2 — sin a ipi +cosa ip2, 
ip' 3 = cos a ip3 ~ sin a ip^ , ip'^ — sin a ip3 + cos a ipi. 

(144) 

This is easily seen from the infinitesimal transformation 

Sip-y = cd-ysips (145) 

and the relations P = -1, I T = -I, [I,T k ] = 0. These 
rotations carry over to the one-particle and one-hole wave 
functions g 7 and q 1 . 

Correspondingly, the complex two component wave func- 
tions ip and x transform as 



tpr = e"tp , X = e "*x- 



(146) 



The 5*0(2) rotations are now realized as U(l) phase rota- 
tions. If we define for a general complex field rj the charge 
Q by the transformation 



n 



e zaQ r) 



(147) 



we infer that ip carries charge Q = 1, while \ nas opposite 
charge — 1. If ip describes degrees of freedom of an elec- 
tron, x describes the corresponding ones for a positron. 
We note that charge eigenstates exist only in connection 
with a complex structure. The real wave function q 1 can 
be encoded both in ip and \ and may therefore describe 
degrees of freedom with opposite charge. 

The action (|30l) is further invariant with respect to dis- 
crete symmetries. Among them a parity type reflection 
maps 



with 



ip(x) -> Vip(x), 



(148) 



(149) 



and we note V 2 = — 1. We may associate V with a CP- 
transformation for Weyl spinors. The action ([30)l is further 
invariant under the discrete transformation ip —> lip. It is 
therefore also invariant under a parity type transformation 
where V is replaced by IV. 



VIII. COMPLEX STRUCTURE 

We have formulated the classical statistical description 
of a quantum field theory for Majorana spinors in terms 
of a real Grassmann algebra. All quantities in the func- 
tional integral (j40|) and the action (|30|) or (l6|)- (fT0|) are real. 
The introduction of complex Weyl spinors in eq. (|1 24[) or 
(|129|) reveals the presence of a complex structure in this 
real formulation. Such complex structures are the basis 
for the importance of phases in quantum mechanics. The 
Schrddinger equation for a one-particle state can be for- 
mulated in terms of a complex wave function, and this 
generalizes to multi-particle states. 

1. Complex Grassmann variables 

Complex Grassmann variables may be introduced in 
analogy to eq. (I124|) 

Ci = -j=(ih + #2) , C2 = -^(ipa + iip4). (150) 
Together with the complex conjugate Grassmann variables 

cr = ^(v-i - #2) , ci - ^(^3 - #4) (151) 

we have for every x and t four independent Grassmann 
variables £i , Ci , C2 , C2 1 which replace tpi,ip2,ip3,ipi- 

A general element of a complex Grassmann algebra can 
be expanded as 



<?c — y ] Cg 1 ...a k ,a 1 ...a l ( 



Xi ...X k ,Xl. ..Xl) 



k,l 



{ a M---U(xk)QM)---Q l (xi)> (152) 

with complex coefficients c. For a real Grassmann algebra 
the coefficients c are restricted by g* = g. From a general 
complex g c we can obtain an element of a real Grassmann 
algebra as 



9 = \ (9c + 9*)- 



(153) 



The action (fBTJf is an element of a real Grassmann algebra 
and reads in terms of £ 

s = [ {c f (d t -T k d k )t+c T (d t -T* k d k )e} 

Jt,x 

= 2 f C\d t -T k d k )C. (154) 

Jt,x 

Up to the factor 2, which may be removed by a rescaling 
of C, this is the action for a free Weyl spinor. 

On every factor £ the action of a matrix multiplication of 
ip by I amounts to a multiplication with i. More precisely, 
we can interpret eq. (|150|) as a map ip — > ^[ip] with the 
property Q[Iip] = iC,[^P\- For the complex conjugate one 
has C*[/V ; ] — ~~K*[4>]- F° r the infinitesimal transformation 
(|145|) one concludes 

<5[C ai (zi)...G ( (zO] =iaQ{ ai (x 1 )...Cc, l (x l ) ) (155) 
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where the charge Q counts the number of factors C minus 
the number of factors for a given term in the expansion 
()152j) . In other words, products with Q + fa_ctors_C and Q- 
factors C* are charge eigenstates with Q = Q+ — Q-. States 
with a given Q are degenerate since many different choices 
of Q+,Q- lead to the same Q. For a real Grassmann al- 
gebra we use eq. (|153[) . Every term in g c with Q =^ is 
accompanied by a term with opposite charge — Q in g*. 
Thus the expansion (|152[) of eq. (|153[) involves "charge 
eigenspaces" with pairs of opposite axial charge. 

It will often be convenient to use an extension to a com- 
plex Grassmann algebra, where the coefficients c in eq. 
(|152[) are arbitrary and charge eigenstates belong to the 
Grassmann algebra also for Q =t 0. This can be mapped 
at the end to a real Grassmann algebra by eq. (|153[) . The 
action of Q± on g c is represented as 



f d 



°-^ c(! "55M' (156) 



The Grassmann derivatives d/d(, a obey the standard anti- 
commutation relations 



d 



d( a (x)' 



>0(y) ? = S a pS{x - y) 



d 



d( a (x) 



(157) 



In the complex basis we can write the time evolution 
equation in the form 



id t g c = Ug c 



with 



H = i 



{TkdkQc 



9C*a 



(r* k d k C) a 



(158) 



(159) 



Both Q+ and Q- commute with %. The operators creating 
one hole or one particle states read 



1 



V2 



d ltd d \ . . 

^W, = VtV a tK« +iPa dC*)- (160) 

The infinitesimal Lorentz transformations of the com- 
plex two-component spinors £, = — ie„ m £ m ™£, are rep- 
resented by the complex 2x2 matrices £ mn given by eq. 
(|127j) . We also observe that £ = EC, = —it-iC, transforms as 



je mrl {2Ji'""" y C such that 



|c r C = ClC2 = \ - lM»4 + i(^A + ^2^3)] (161) 

is a Lorentz scalar. The same holds for QQ such that 
the bilinears rpi?p3 — V'2'04 and ipiipi + ip2 "03 are separately 
Lorentz scalars. Also the product 

C1C2C1C2* - M2^4 = ^e a ^ s ^^ s (162) 



is a Lorentz scalar. The functional measure obeys 

dip4dip3dip2dipi = ]^[ dx/j-y 
1 

J d( 2 dCi J dQdCx = J dCd(*. (163) 



2. General complex structure 

A real Grassmann algebra will, in general, admit differ- 
ent possible complex structures. A discussion of phases in 
quantum mechanics needs a specific choice of the complex 
structure. For example, our discussion of Dirac-spinors in 
sect. |X]will employ a complex structure different from eq. 
(|150p . We therefore briefly discuss the general properties 
of complex structures. 

In a real even-dimensional vector space a complex struc- 
ture is given by the existence of an involution K, together 
with a map /, obeying 



K 2 = 1 , I 2 = -1 , {K,I} = 0. 



(164) 



The matrix K has eigenvalues ±1 and we may denote eigen- 
states with positive eigenvalues by vr and those with neg- 
ative ones by vi,Kvr = vr,Kvj = —vj. The matrix / is 
a map between vr and vi, implying that the number of 
independent vr and vj are equal. We can choose the vj 
such that Ivr — vj, Ivj — —vr and use these properties for 
defining a map from the real vectors v to complex vectors 
c = vr + iv 1 = c(v), with the properties c(Kv) = [c(v)\ , 
c(Iv) = ic(v). A linear operator or observable A can 
be represented by multiplication with a complex matrix 
if [A, J] = 0. 

As an example, we consider the complex structure under- 
lying eq. (IT50| . The maps (fT24|) or (IT50| act onu = {g 7 } 
or v = {ip-y } , and the matrices K, I are given by 



K=K=\ ° 1 . T = I. 



(165) 



Indeed, the map ([tp] given by eq. (I150|) obeys 

c[^] = rw,cN=iCM. (166) 

A transformation ip — > Aip is compatible with this com- 
plex structure if A obeys 



[A,I}=0. 



(167) 



In this case A can be represented by complex matrix mul- 
tiplication acting on 



CL40] = IC[# 



(168) 



The matrices T k commute with / and are therefore com- 
patible with the complex structure 



CPW] = t*C[# 



(169) 



Also / and (cf. eq. (|3"9"|0 are compatible with the 
complex structure, where the action of E M1/ on £ is given 
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by eq. (|127l) . The matrices A obeying eq. (jl6T[) form 
a group, since the product of two matrices AiA 2 again 
commutes with /. This product is represented by complex 
matrix multiplication of A\ and A2, ^[AiA^ip] — A\A2C[ip]- 
In contrast, the matrices 7 M anticommute with /, cf. eq. 
()3GI) . and are therefore not compatible with the complex 
structure (|165[) . One finds 



C[ 7 >] = -T2C 



(170) 



The complex structure (|164[) , (II 65[) for the Grassmann vari- 
ables can be extended to a complex structure for the real 
Grassmann algebra by defining a suitable map g — > g c . 



IX. MASSIVE MAJORANA FERMIONS 

The action ([50)1 describes massless Majorana or Weyl 
spinors. For a massive Majorana spinor one adds a mass 
term 



ip(ml — m)ip 



^ 7 [m( 7 /) 75 -m( 7 °) 75 ]^. (171) 



In the presence of the mass term (|171|) the action remains 
Lorentz- invariant, anti-hermitean and real (for real m, to), 
and the Minkowski action Sm = iS is hermitean. For 
discrete time steps the first factor ip is taken at t, and the 
second at t + e. The E/(l)-rotations (|144p do not leave S m 
invariant. Indeed, the infinitesimal transformation (|145[) 
results in 

M -> M + SM , SM = -a[I, M] = 2aMI, (172) 
where 

M y5 = m(j°I)^s - to( 7 %-. (173) 
The terms ~ m and to are rotated into each other, 



Sm = —2am , Sm = 2am. 



(174) 



We may use this transformation in order to choose a con- 
vention where m = 0. With respect to the parity transfor- 
mation V (|149|) the term ~ m is invariant, while m changes 
sign. On the other hand, S m is odd with respect to the 
discrete transformation ip — > lip. Thus, with respect to 
the parity transformation IV one finds that m is invariant 
while to changes sign. 

1. Evolution equation 

For massive Majorana spinors the Grassmann time evo- 
lution ([57)1 obeys 



K 



with 



d 



x dip 7 (x) 



{(Tk)jsdk~ M yS }ips(x) =/C +/C m , (175) 



d 



dlpy(x) 



M jS ip s (x). 



(176) 



The particle number Af remains conserved, while the 
charge Q is no longer a conserved quantity in the presence 
of a mass term. 

We may consider general operators of the type 



B ev (y) 







Be V (y)ipr,{y), 



(177) 



with B t1) depending on y and derivatives with respect to y. 
They obey the commutation relation 



[Ben(y),t£o] = d k |^-(T fc ) 7£ B er ,V t) 



(178) 



d 

+ -Q^-{B lt {T k ) eil dk^ v - (T k ) 7> :dk{B en ip v )}, 

where all quantities on the r.h.s. depend on yk and dk = 
d/dyk- For the mass contribution one finds 



[B en (y),K, m ] = 



d 



dtpj(y) 



[M,B(y)] jS My) 



(179) 



The momentum operator (II 26[) reads in the real Grass- 
mann algebra 



f d 

Jy Ol/J 



(180) 



For M 7^ it is no longer conserved, while we may define 
conserved quantities involving an even number of deriva- 
tives 



d_ 

y dtpj 



(181) 



For any static vacuum state ga,)Cgo = 0, the evolution 
of the one-particle state (|116[) follows eq. (|120j) 



d t q 1 {x) 



d 



dip 1 {x) 



9o 



q 1 {x) 



8 



dip 7 (x) 



go- (182) 



With 



d 



8ljj-y(x) 



9 »V 9 

-M 5l =M lS - 



d^ S {x) 07 70 dll> S (X) 

one finds the evolution equation 

9tQ = (T k d k - M)q , Vq = ( 7 ^ M - M)q = 0, 
where 

M = - 7 °M = ml — to. 
For to = the relation { 7 M , 1} = implies 



By using 



V 2 = 9^ -to 2 . 



[K m ,ip-y(x)] = M^ s ip s (x) 



(183) 



(184) 



(185) 



(186) 



(187) 
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we find for the one-hole wave function q (111 8[) the same 
evolution equation as for the one-particle wave function q, 



Vq = 0. 



(188) 



It is instructive to write the one-particle evolution equa- 
tion (|184l) in terms of the two-component complex Weyl 
spinor ip defined by eq. (|124[) 



d t ip = T k d k <p - iT 2 (m - irh)p*. 



(189) 



The fact that ip* appears shows that the multiplication 
with M is not compatible with the complex structure (|165|) . 
In terms of \ (|1 29[) and ipM (|131l) one finds 



d t p = T k d k (p + (m - im)x, 
dtx = ~T k d k x - (m + im)tp , 



(190) 



and 



d t ipM = 7(V)(7(W) 5 fe ~ im l(w) + ™)ipM, 
(itfwfa + m l(w) + im)ip M = 0, (191) 

with Trejn the Dirac matrices in the Weyl representation 

urn 

2. General solution for single massive Majorana 
particle 

In order to find the general solution of eq. (|184[) we 
expand similar to eq. (|13T[) . We may consider modes with 
a fixed momentum p. (As familiar in quantum mechanics, 
these modes are not normalizable for infinite volume and 
may be considered as limiting cases of normalizable wave 
packets.) For every given p ^ we can define 



H(p) = 



\p\ 



H 2 (p) = 1 , PkPk = \p\ 2 , 



(192) 



such that H(p) has eigenvalues ±1. Decomposing (p{p) 
according to the eigenvalues of H(p), 



and using 



H(p)p±{p) = TV±{p), (193) 



X±(t,p) = -iT2tp±(t,-p), 
H(p)x±(t,p) = TX±(t,p), (194) 



one has 



p{t, x ) = / {(p+^py^ + cp^py**), 
jp 

X (t,x) = f (x+(t,p)e tpx + X-(t, P y px ). (195) 



This yields, for to = 0, 



dt<P±(p) = Ti\p\<f±(p) + mx±(p), 

dtX±(p) = ±i\p\x±(p) ~ m<P±(p), (196) 



with general solution 

<p+(t,p) = 6 + (p)e-' wt -^T 2 6;(-p)e« wt , 

to 

X+(t,p) = -ir 2 b* + (-p)e iut - i^Mb+We-™ 

to 

<p-{t,p) = b-ipy^ + ^^T 2 b*_(-p)e-^ \ (197) 

m 



X-{t,p) = -ir 2 b*_(-p)e- wt + i ^b.{p)e^\ 

to 



where 



= w(p) = +v>l 2 + ™ 2 - 



(198) 



3. Classical probabilities for Majorana fermions 

The general solution for tp(t, x) depends on four complex 
functions contained in b± (p) . From the real and imaginary 
part of <p(t, x) we can infer the four real functions g 7 (x) and 
therefore the classical probability distributions which de- 
scribe the one-particle states. As an example, we consider 
the special solution b-{jp) — 0, b + ^(p) = 0,& +i i(p) = f(p) 
real, for which 



r - j( /.,i = / ^^Mfi-py^), 

to 



(199) 



qi(t,x) = / f(p)cos(px-u>t), 
Jp 

q 2 (t, x) = I f(p) sin{px - ut), 
Jp 



Qa(t,x) 



f{p)- ^- sin(px - ut), 

to 

f(P) U ~ * cosOpx - ut). (200) 

TO 



We may further take f(p) = fS(p — p) with p = 
(0, 0, -P),P > 0, such that, with u) = y/P 2 +m 2 , 

qi(t,x) = / cos(Px 3 + tot) , q 2 (t,x) = -f sm(Px 3 + Lot), 

LO — P 

qa(t,x) = / sin(Px 3 + tdt), 

TO 

LO — P 

qi{t,x) = -f cos(Px 3 +ujt). (201) 

TO 

For this particular example the normalization J" q 2 (x) = 1 
requires (with V the space-volume) 



1 



2V m 2 + P 2 - PVP 2 + m 2 



(202) 



For the general solution (|200l) the probabilities p 1 (t, x) = 
q 2 (t,x) describe the time dependent probability distribu- 
tion of a classical statistical ensemble which accounts for 
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all properties of a propagating massive single Majorana 
fermion. If the vacuum is the totally empty state g$ = |0) 
the states (7, x) correspond to "Ising states" r where only 
one bit at x of species 7 is occupied and all others are 
empty. (For a different go the classical states associated 
to (7, x) are more complicated.) The probability of a state 
(7, x) is given by p-y(x) = q^(x), e.g. according to eq. (I200p 
with a proper normalization of f(p). 

The probability distribution |p 7 (a;)} properly accounts 
for all characteristic quantum interference phenomena for 
modes with different momenta p. This extends to the 
correct description of interference for mult i- particle states. 
(For a concrete example of a two-particle state for massless 
Majorana spinors cf. ref. ||.) The appearance of interfer- 
ence phenomena is not surprising in view of the formulation 
of the time evolution in terms of a Grassmann functional 
integral. On a deeper level, it is related to the fact that sim- 
ple time evolution equations ([3]) are linear in the classical 
wave function q T , but not linear in the classical probabili- 
ties p T . 



X. DIRAC FERMIONS 



express the action fl203[) in terms of if> o as 



where we define 



ipD^df, - mI)ip D , 



4>d = 0^7°- 



(206) 



(207) 



The matrices 7 M and / appearing in eqs. (1206[) . (|207|) are 
all real. We may also introduce a purely imaginary matrix 



*7°7 1 7 2 7 3 = —H j T 4 = 1. 



(208) 



It is hermitean and anticommutes with all Dirac matrices 



7 f =7, {7,7 M } = 0, [7, EH 



0. 



(209) 



Often 7 is denoted as 7 s . On the level of the eight variables 
0i i7 and "02,7 the matrix 7 is defined as a real 8x8 matrix, 
with I2 not acting on the index 7, 



Dirac spinors can be composed of two different Majorana 
spinors with equal mass. The action reads 



S = 



{Vi(5 t -T fc 5 fc + m 7 °/> 1 



+^2(ft-r fc a fc +TO7 /)V 2 }, 



(203) 



with Tfc, /, 7 given by eqs. (JTT1), ©, ([3"I]). (For simplicity 
we take in — in eq. (HI]).) The functional integral 
extends now over the Grassmann variables 0i l7 (i, x) and 
ip2,j(x,t). For a free theory the action involves separate 
pieces for ipi and 02) S = Si + 82- If the initial state 
g(to) = gi n factorizes into two factors, one involving only 
01 and the other only "02, the Grassmann wave function 
g(t) factorizes for all t. There is no need, however, for 
g(to) to factorize and general states g(t) will not be of a 
factorisable form. 



7 = —h ® I i h = 



0, -1 

1, 



(210) 



(The matrix / defining the complex structure is given by 
I2 according to eq. (|205p .) 

Since 7 commutes with the Lorentz-generators E Ml< the 
projections 



01 



1 



(1 + j)lpD , Ipi 



1 



(1-7)^ 



D 



(211) 



are representations of the Lorentz group. They transform 
as the two inequivalent fundamental spinor representations 
and describe Weyl spinors. In terms of ipi,>y, ^2,7 one has 



L.R 



( "01,1 T ^2,2 + #2,1 ± #1,2 

"01,2 ± ^2,1 + #2,2 T #1,1 
01,3 T 02,4 + #2,3 ± #1,4 
V 0)1,4 ± ^2,3 + #2,4 T #1,3 



(212) 



1. Complex structure 

We introduce a complex structure by defining the four- 
component complex Grassmann variables tpu 



■0D = "01 + #2 , ipn = ipi- #2 



(204) 



In terms of the general discussion in sect. IVIIII the matrices 
K and I act as 



K 



V>l\ #01 

-1P2J ' V02 



"02 
0W' 



(205) 



such that 0i is the "real part" of 0d, and 02 the "imag- 
inary part" . The complex structure employed for Dirac 
spinors differs from the complex structure (|165[) . We can 



where the upper sign relates to 0l and the lower to ipR. 
We observe that both tpL and 0_r receive contributions both 
from ipi and 02- They differ from the Weyl spinors CiiC2 
that can be formed from ip\ and 02 according to eq. (|150[) . 
Using 7 we can write the action (|206[) as 



Sm — — 



0d(7' 1 ^ - imy)ipi 
0z)(#' i d M + W7)0z 



M ■ 



(213) 



Defining the hermitean conjugation as a map ipo ~ > 0|j, 
accompanied by a complex conjugation of all coefficients 
in the Grassmann algebra and a total transposition (re- 
ordering) of all Grassmann variables, the hermiticity of the 
Minkowski action, S\,j — Sm, is easily verified. 
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2. Charge and electromagnetic fields 



The action (|203|) is invariant under all symmetry trans- 
formations defined in the preceding sections, i.e. those 
leaving S\ or S% separately invariant. There are further 
symmetries related to transformations between ipi and tp2- 
Rotations between ipi and ip2 can be accounted for by the 
infinitesimal transformation 



dip 



= -Phi> = 



(214) 



The discretization is the same for ipi and tp2- It therefore 
respects the invariance with respect to the transformation 
(|214j) . This rotation can be transferred to a phase rotation 
of the complex Dirac spinor 



SipD = -i(3*l>D- 



(215) 



We associate the transformations (|214[) . (I215|) with the 
[/(l)-transformations of electromagnetism. The Dirac 
spinor carries charge Q = — 1 and can be identified with 
the electron. Both components ipL and ipn (|211l) transform 
as Weyl spinors with the same electric charge Q. 

We may extend the action by considering Dirac spinors 
in an external electromagnetic field A^{t, x). As usual this 
is done by replacing the derivatives in eq. (|213[) by 



covariant derivatives 
action a piece 



ieA u . This adds to the 



AS = -ie / ipD^A^o. 



(21( 



In terms of the real Grassmann algebra involving ip\ and 
ip2 the additional piece reads 



AS = -e {^(Ao - A k T k )^ 2 ~ M A o - A k T k )^}. 

Jt,x 

(217) 

For discrete time steps the first Grassmann variable is 
taken at t and the second at t + e. In the continuum limit 
the two terms in eq. (|217l) are equal. We may consider 
the fields A^ as "sources" for fermion bilinears involving 
ipi and tp2 ■ For nonzero A^ the action is no longer a sum 
of independent pieces Si + S2 for ipi and ip2, respectively. 



3. Weyl representation 

One may use the similarity transformation (|135[) 



ipD,w = A , 7 ( A 



w) - A 



^(M) ' 



(218) 



in order to bring the Dirac matrices to the Weyl repre- 
sentation (|132[) . The action (|213[) . (|216l) retains its form, 
with the replacements 7^ — > Jn^y 4>d — > 4'd,w,''Pd 

= 4>d,w- In the Weyl representation the ma- 
trix 7 is diagonal, 7 = diag(l, 1,-1,-1). The Dirac spinor 
tpD w can therefore be written in terms of two-component 



spinors ipL and tpR, 

WD,W 
Ca 



, tpL = (Cl + K2) , ^R = (Cl + 



^R 

1 flpa, 1 + ilpa,2 
V2 ^0,3 + ilpa,iJ ' 



t 1 f-1pa,3 + I1pa,4\ . A * 

= —7={ , ■ 1 ) = - lT ^a> 

V2\ Va,l-Wa,%J 



(219) 



which, in turn, can be written as linear combinations of Cl 
and C2 or Ci and £2- 

We observe that in this basis the operation of complex 
conjugation is represented by an involution K different 
from K in eq. (|205[) . In terms of the eight-component field 
ip = (ipi,4>2) one has K = diag(l, — 1, 1, — 1, — 1, 1, — 1, 1). 
The matrix / remains I2 as in eq. (|205|) . and I2 and K 
anticommute 



K 



K, 
0, —K 



0, -1 

1, 



{K,I 2 } = 0. (220) 



Thus the similarity transformation (I218P changes the com- 
plex structure. On the level of a real representation it can 
be expressed by real 8x8 matrices A' such that a similar- 
ity transformation may, in principle, transform both K and 
/. The matrix / remains unchanged, since iAtp — A(itp) 
in the complex representation. For A ^ A*, however, the 
matrix K changes since (Atp)* ^ Aip* . 

In the Weyl basis (|219[) we can take ipL and ipR as two- 
component complex spinors. In terms of ipL and i/jr the 
action reads 

S = - {^L^id^ + ieA^L+rPnYid^ + ieA^jR 

Jt,x 

-im^RipL-ipM}, (221) 
where we define 

$ = {^R,^l) = -i(ip R , V>1) 
4>l = ^ ,t>R = t> (^) ■ (222) 

For m = it is invariant under separate "chiral phase 
rotations" of ipL and ipR. The transformation 5^>l = 
iaipL^ipR = —iaipR corresponds to 5( a = iaS( a and there- 
fore to a simultaneous transformation of the type (1145[) for 
ipi and ip2- 

4. Left-handed representation 

We can also reformulate the action in terms of two left 
handed Weyl spinors ipL and ipi by introducing 

Ipt = lT2lp R = Cl - <2, 

iPr = -iT 2 (r L )* , 4>R = {^ c l) T T2. (223) 

This yields 

S = WlI 3 * + ieA - T k (d k + ieA k )]ipL 
Jt,x *• 

+(^£)t [8 t - ieA Q - T k (d k - ieA k )]i,l 



(224) 
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The charge of ipl is opposite to tpL - if ipL describes a left- 
handed electron, -0£ accounts for the left-handed positron. 

We may group ipL and ipl into a four-component com- 
plex spinor 



LL 





(Ci- 


MC 2 \ 









(225) 



In this representation the matrix / defined by eq. (fl3|) is 
represented by a multiplication with i, 



1pLL[Ilp] = ilpLL, 



(226) 



while the matrix I2, which acts on (Ci ? C2) or (V^V^) as 
defined in eq. (|210p . is now represented as 



4>ll [hip] 



i, 
0,-i 



LL- 



(227) 



We observe that the role of I and I2 is exchanged if we 
compare the complex spinors %/jll and ipD,w> 



1pD,w[hi>] = ilpD,W, 

'i, 



(228) 



ipD,w[Iip] 



0. 



i>D., 



W 



D,W- 



The complex structures associated to ipLL and ipD,w are 
therefore different. Nevertheless, the operation of com- 
plex conjugation of ipLL or ipD,w corresponds in the real 
representation ip to the same matrix multiplication with 
K = diag(K,~K). 

5. Enhanced symmetries for free massless Dirac 
spinors 

For m — 0, — the free theory exhibits a symmetry 
of S't r (2)-rotations among ipL and ^>£, 



8tp 



LL 



iUjTjljjLL, 



(229) 



with 0J3 = j3 identified with the transformation (I215|) and 
Tj the Pauli matrices acting on the two components of ip2d 
(not on the spinor components of ipL and For exam- 
ple, the transformation with u)\ acts as 

SlpL = -UJlT 2 lp* R , SlpR = -UJlT 2 1p* L , 

SCi = iwiCi , SC2 = -iwiCj. (230) 



This corresponds to a transformation of the type (|145|) . 
now with opposite directions for ip\ and ip2 ■ The situation 
is similar for the transformation ~ u>2, where 

Sd = -iui 2 (2 , K2 = -«w 2 Ci, 

Sipi = ~u} 2 Iip2 , Sip2 = (231) 

leads again to a transformation among ipi.-y and -02.7 which 
involves /. 



6. Evolution equation for Dirac fermions 



For the action S + AS f2"05|) . (f2T7) the Grassmann evo- 
lution equation dtg = /Cg involves K = /Co + /C m + AIC 



/Co + /C m 
AK, 



E 



iTA-mj !)^), 



d 



d 



(A (x) - A k (x)T k )Mx) 



dip 2 {x) 



(A (x) - A k (x)T k )Mx) 



.(232) 



Here we work in the Majorana basis (fTTj) . (fT3|) . (|31|) and we 
have suppressed the spinor index 7. The evolution equation 



d t g(t) = Kg(t) = (/C + K m + AK)g(t) 



(233) 



describes the dynamics for an arbitrary number of charged 
relativistic fermions (and their antiparticles) in external 
electromagnetic fields. 

The evolution equation for the classical wave function 
(1861) . and therefore for the probability distribution, ob- 
tains from eq. In summary, we have formulated a 
time evolution equation for a classical wave function q T (t) 
which describes an arbitrary number of charged electrons 
in external electric and magnetic fields. In the next section 
we will see that its restriction to one-particle states yields 
the relativistic Dirac equation. As usual, a non-relativistic 
approximation will yield the familiar Schrodinger equation 
for a particle in a potential or moving in external magnetic 
fields. 

7. Classical observables 

Classical observables can be constructed from linear 
combinations of products of occupation numbers. We com- 
bine the index a = 1,2 for the two Majorana spinors ip a 
with the index 7 into a common index e = (7, a), e = 1 . . . 8. 
The Grassmann operators corresponding to the occupation 
numbers N e (x) read then 



d 



(234) 



dA(x) 

In a formulation with discrete lattice points they obey 

M?{x) = K{x) , [K(x),M n (y)] = 0. (235) 
The basis elements g r are eigenstates of Af e (x), 

N,(x)g T = (N e (x)) T g T , (236) 

with (N e (x)J = if g T contains a factor ip e (x), and 
(iV e (x)) = 1 if no such factor is present. The expecta- 
tion value of a product of occupation numbers obeys 

(N €l ( Xl )N e2 {x2)...N em (x m )) 

V<ij){t)g(t)N ex {xx) • . .K m (x m )g(t) 



Y / Pr(t)(N ei (x 1 )) T . . . (N em (x m )) T , (237) 
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corresponding to the standard rule for classical statistics. 
Obviously, the expectation values of classical observables 
can be computed from the probability distribution p T (t) 
and do not involve the signs s T in eq. ©. Classical prod- 
ucts of classical observables commute. 

In order to compute the time evolution of expectation 
values of classical observables we define 



M e n(x) 







dip e (x 



(238) 



and use the relation 



[M ev {x),K] = X)tt T *)«A 



■ip v (x) 



d 



-d k tp a (x) 



Jv "dMx) 
+W ea (x)M av (x) - M ea (x)W av (x)} 



(239) 



with 



W ev (x) = -m(~f° I) al 35 ab + e(A (x)S al 3 - A k (x)(T k ) a/3 )e ab , 

(240) 

and d k = d/dx kl e ab = -e ba ,e 12 = l,r) = (/3, 6), (T k ) ev = 
(T k ) a p5 ab = (Tfc)^ e . We observe that W is antisymmetric, 
W er] (x) = — W ne (x). A local particle number can be defined 
as 



M(x) = j2^(x) = J2 M ^ x )> 



and obeys 



[Af(x),K] = {T k \ a d k M 



cx q • 



(241) 



(242) 



According to eq. (|107l) the time evolution of the classical 
mean local particle number involves expectation values of 
"off-diagonal" operators M mi 



The particle number 



N = / Af(x) , [Sf, fC] = 0, 



(243) 
(244) 



is conserved. 



XI. QUANTUM MECHANICS FOR PARTICLE IN 
A POTENTIAL 

1. Dirac equation 

We next concentrate on one-particle states, described by 
Grassmann elements 

9i(t) = [ ( Qi,j(t, x) - + g 2 , 7 ft, x) 



<9i/>i, 7 (x) 



dip 2 ,~,(x) 



.90- 

(245) 

Here go is some arbitrary static "vacuum state". For 
^-So = we find for the one-particle wave function the 
Dirac equation 

(7^ - ml)qi = eA^q 2 , 

(Vfy - ml)q 2 = -eAyrfqx . (246) 



For the derivation of the Dirac equation we use 
d 



AK. 



dtpi n (x) 
e ( Ao(x) 



d 



9^2,7 (x) 



A k (x){T k ) Sl 



d 



dil)2,s{x) 



A/C 



d 



dip2,~f{x 
e I A Q {x) 



(247) 







dipi n {x) 



A k (x)(T k ) Sj 



d 



dtpi,s(x) 



The time evolution equation for the one particle wave func- 
tion then extends eq. f|l 20[) to 

d t qi = {T k d k - m 7 °/) gi + e(A - A k T k )q 2 , 

d t q 2 = (T k d k - mi I)q 2 -e{A -A k T k ) qi . (248) 



This is equivalent to the Dirac equation (|246l) which can 
also be written as a matrix equation. 

( 7 ^-mJ)Q =0, 

D M = d^ + eA„ Q = d» + eA^h . (249) 

The usual complex form of the Dirac equation is recovered 
if we use a complex one-particle wave function 

<Pd — qi + iq2 , 7 M (^ + ieA^tpD = im^ip D . (250) 

This equation holds for an arbitrary representation of the 
Dirac matrices 7^. 

The derivation of the Dirac equation for a one-particle 
state (|245[) has only used the general evolution equation 
for Grassmann elements g(t) and the condition JCgo = 0. It 
therefore describes the dynamics for a very extended fam- 
ily of classical probability distributions or classical wave 
functions. Indeed, it is sufficient that go is an arbitrary 
static state (not necessarily a priori with a fixed particle 
number). This reflects the physical property that isolated 
one-particle states can occur under a wide variety of cir- 
cumstances, and that for sufficient isolation the proper- 
ties of the environment do not matter for the dynamics of 
the isolated particle. In this context we emphasize, how- 
ever, that our model only describes external electromag- 
netic fields while we do not account for the fields generated 
by the particles that may be present in the environment. 
In this sense our setting describes "real physics" only in a 
situation where the electromagnetic fields generated by all 
present particles can be approximated by a "mean field" 
that is independent of individual particle positions. This 
is precisely the setting of standard one-particle quantum 
mechanics. 

The presence of electromagnetic fields influences the con- 
ditions for a static "vacuum" or "environment", ICgo = 0. 
We will concentrate here on simple vacuum states as go = 1 
or go = |0), where |0) involves a product of all Grassmann 
variables and obeys Af\0} = 0. These states are static for 
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arbitrary electromagnetic fields. For simplicity we choose 
9 a = |0) such that 



M 91 



9i- 



(251) 



(The following discussion can be easily generalized, how- 
ever, to all static states go which are eigenvalues of 
Af,Afgo = mj , by using a modified local occupation num- 
ber where a suitable constant is subtracted from Af (x).) 

2. Particle observables 

Let us consider the expectation value of the local occu- 
pation number 



We use the identity 



8 



dip* (x) 



in order to establish the simple relation 
(A^(x))=]>> 2 (i,x). 



(252) 



(253) 



(254) 



This holds for a wave function g(t) = gi{t), and we use 
eq. (|103|) employing appropriate basis elements (a subset 

0i{9r}) 



9iA x ) 



d 



h,ri(y) = i>ri{y), 



dip e (x) 



(255) 



with 



9i = / qe(x,i)g li€ (x) , gi= q v {y,t)gi, v {y)- (256) 

J x J y 

The normalization J Vipgigi = 1 implies for the one- 
particle wave function the normalization 



J X J X 



(257) 



which is compatible with (M) = 1, cf. eq. (|251|) . The nor- 
malization (|257l) is preserved by the unitary time evolution 
of the one-particle wave function. Indeed, we can write the 
Dirac equation (|250p as a Schrodinger-type equation 

ihdtfD = Hip D , (258) 

with hermitean Hamiltonian (Tk = j ^) 

H = ihT k d k + ftm 7 °7 + he(A - T k A k ) - if. (259) 

For a pure one-particle state the expectation value 
(N(x)} can be interpreted in a natural way as the prob- 
ability density to find the particle at the position x. This 



is precisely the standard interpretation of i/ , £)(x)V'd (x) in 
one-particle quantum mechanics, 



w ( x ) = (N(x)) = ip f D {x)ip D (x), 
w(x) = 1. 



(260) 



Since N(x) is a classical observable, we can define the po- 
sition of the particle as a classical observable 



X = / xN(x). 



(261) 



The expectation value in classical statistics coincides with 
the standard quantum mechanics rule 

(X) = ( / xN(x)) = I xw(x) = I (p]j(x)x(p D (x). (262) 



This finds a natural extension to arbitrary functions of 
the position observable 



f(X) = / f(x)N(x). 



(263) 



The expectation values of this type of classical observables 
follow again the rule of quantum mechanics 



(f(X)) = / <pUx)f{z)<p D (x). 



(264) 



In particular, one obtains the same formula for the disper- 
sion (XfcXfe) — (Xk)(Xk) as in quantum mechanics. We 
conclude that measurements of the position of a particle, 
or more generally the distribution of positions in an ensem- 
ble of one-particle states, can be described equivalently in a 
classical statistical ensemble with classical observables, or 
in quantum mechanics with Hamiltonian (|259j) . The com- 
plete time evolution of the distribution of positions is iden- 
tical in both descriptions. This covers, in particular, the 
characteristic quantum interference in a double slit experi- 
ment. The time evolution (|86p . (|52"|) . (|233[) for the classical 
wave function and associated classical probability distri- 
bution {p T } produces for one-particle states exactly the 
quantum mechanical interference pattern. 

3. Schrodinger equation 

Standard quantum mechanics for an electron in a poten- 
tial is recovered from the non-relativistic approximation to 
the Dirac equation. This is well known, and we sketch 
here for completeness only the case A k = 0. The non- 
relativistic approximation becomes valid if eAo and iT^dk 
are small compared to m. Since (7°7) 2 = 1, it is con- 
venient to choose a basis where 7°7 = diag{\, 1,-1,-1). 
With V(x) = HeA Q (x) and M = hm the Hamiltonian (|259l) 
takes the form 



M , a kPk \ 
H = I t +V{x), 
a kPk , -M J 



(265) 



where we use the standard quantum mechanical momen- 
tum operator 



p k = -iWk- 



(266) 
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The matrices 

(Ti = -1 , a 2 = ir 2 , <t 3 = ir 3 (267) 

arise from 

T> k = tfT k U = -^ > (268) 

with 

u = ±(X ; ,c/t C7 = 1)C/ t = C ;. (269) 

(Note T-f = Ti, = -T 2 .) The unitary matrix U diago- 
nalizes 7°7, 

We observe for the matrices at the relations 

a \ a i + a \ a k = o k a\ + dial = 24;, (271) 

which guarantee {I^T/} = 28u- 

We next decompose tpo into two-component wave func- 
tions, ip'jj — [x T ,P T ), which obey 

ihdtx = (M + V)x + <JkPkP, 

ihd t p = {-M + V)p + 4 PkX . (272) 

For the non-relativistic electron we consider the approxi- 
mate solution p = A\, where A = a\pk/(2M) is deter- 
mined by requiring in leading order dtp = AdtX- Insertion 
into eq. ([21]) yields for tp = exp(iMt/h)x the standard 
Schrodinger equation for a particle in a potential V, 

W=g+V)^. (273) 

As usual, the Dirac equation can also describe non- 
relativistic positrons. For this purpose one considers a sec- 
ond class of solutions x = Bp, with B = —akPk/ For 
ip = exp(iMt/h)p* one obtains 

lhd ^=&^- V )^ (2?4) 

and we note the change of sign of the potential due to the 
opposite charge of the positron. 

All quantum mechanical phenomena extracted from so- 
lutions of the Schrodinger equation are described by our 
time evolution equation for a classical statistical ensem- 
ble of Ising-spins on a lattice. For a potential realizing a 
double-slit situation the standard interference pattern is 
predicted for this classical statistical ensemble to appear 
behind the slits. This holds provided that the initial state 
at some time to corresponds to a one-particle state describ- 
ing a particle moving towards the slits. 



4. Particle-wave duality 

The discreteness of measurement values in quantum me- 
chanics can be traced back to the discrete occupation num- 
bers of the Ising-type model. In quantum mechanics we 
may define an "interval observable" Jn by a function 

J ^)-{Jo f thefwTse)- ^ 
It has the property 

4 = J n , (276) 

such that its spectrum consists of the discrete values and 
1. According to the rules of quantum mechanics, the pos- 
sible outcomes of a measurement of J-r are or 1. The 
interpretation in quantum mechanics is simple: either the 
particle is within the region (interval) J-r, in which case the 
measurement value J-r = 1 will be found, or it is outside 
this region, and Jn = will be found. Particles are discrete 
objects - they are either inside or outside an interval. 

In our classical statistical Ising-type setting J-ji is a clas- 
sical observable, given by 

Jn= I N(x) = Y^N L (x), (277) 
Jn n 

The sum ^ K extends over all lattice points within the re- 
gion 1Z. The classical observable Nl(x) corresponds to a 
normalization of occupation numbers for the discrete lat- 
tice where (Nl,c) — 0, 1. The sum over species Nl(x) = 
J2^L,e(x), cf. eq. (|252[) . can therefore take in any classi- 

€ 

cal state r only the discrete values [Nl(x)) t — (0, 1 . . . , 8), 
according to the total number of eight species. In conse- 
quence, for any state r of the classical statistical ensemble, 
(Jk.)t is a positive integer or zero. According to the stan- 
dard rule of classical statistics these integers describe the 
possible outcomes of measurements. 

For a one particle state the total particle number equals 
one, 

1= f N{x)=Y,Nl{x), (278) 
v v 

where ^2 extends now over all lattice points in the total 
v 

volume. Since 1Z must be contained in V the maximal 
allowed value for J-r in a one-particle state is one, such that 
(Jtl)t = 0, 1 are the only possible values of the classical 
observable for such a state. This is precisely the quantum 
mechanical rule. No new postulate is necessary for this 
measurement in quantum mechanics - the quantum rule is 
inferred from the standard rule of classical statistics. Of 
course, the expectation value of J-r for a one particle state 
is the same in the quantum mechanical and the classical 
statistical description 

(Jn) = / <Pd( x )Jk(x)vd(x) = / (p' D (x)(p D (x). (279) 
Jx Jn 
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While our model connects the discrete particle aspects 
in quantum mechanics directly to the discrete classical 
Ising-spins, the continuous wave aspects also arise in a 
natural way. The quantum wave function is continuous 
because probability distributions and associated classical 
wave functions are continuous (at least piecewise). As we 
have mentioned already, the characteristic interference ef- 
fects for waves and the superposition arise from the linear- 
ity of the fundamental evolution equation in the classical 
wave function. 



XII. CONCLUSIONS 

We have derived the Schrodinger equation for a quan- 
tum particle in a potential from a classical statistical en- 
semble for Ising-spins. The dynamics of the classical sta- 
tistical system has to be specified by an appropriate evo- 
lution equation for the probability distribution. For this 
purpose we employ the classical wave function which is 
defined as the positive or negative root of the probability 
distribution. The proposed evolution equation is a linear 
differential equation for the classical wave function. The 
wave function at time t obtains from the wave function at 
t' by a rotation - this guarantees the preservation of the 
normalization of the probability distribution. 

We have exploited a map between the classical wave 
function and a Grassmann wave function which is an el- 
ement of a real Grassmann algebra. In turn, the time evo- 
lution of the Grassmann wave function can be associated to 
a Grassmann functional integral. This allows us to formu- 
late the evolution equation for the classical wave function 
in terms of the action of a functional integral. The symme- 
tries of the model, as Lorentz symmetry and electromag- 
netic gauge symmetry for our model of Dirac spinors in an 
external electromagnetic field, can be easily implemented 
in this way. Since the classical wave function is real we 
have to formulate the model in terms of a real Grassmann 
algebra. 

Our model is regularized on a lattice of space points. 
On the one hand, this guarantees that mathematical ex- 
pressions are well defined for a finite number of lattice 
points, with continuum limit of an infinite number of points 
taken at the end. On the other hand, the concept of classi- 
cal Ising-spins or associated occupation numbers at every 
point x is well defined. The discreteness of the particle 
aspects of quantum mechanics can be traced back to the 
discrete occupation numbers that can only take the val- 
ues zero or one. We also have used discrete time steps, 
and we have formulated the lattice action such that the 
Grassmann wave function g(t + 2e) obtains from g(t) by 
a rotation. The time evolution is unitary not only in the 
limit e — > 0, but also for finite e. A unitary time evolution 
for infinitesimal time steps is easily achieved by any anti- 
symmetric evolution generator K Tp for the wave function 
Qr{t), 
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Generating a unitary evolution also for finite smallest time 
steps imposes restrictions on the form of the lattice action. 
This, together with the requirement of a real Grassmann 
action, requires some care for the construction of the model 
and explains the specific form of the action in comparison 
with other possible lattice actions. 

The proposed evolution equation for the classical sta- 
tistical ensemble of Ising-spins does not only lead to 
the Schrodinger equation for non-relativistic one-particle 
states. It entails the full dynamical equations for a quan- 
tum field theory of Dirac fermions in an external electro- 
magnetic field. The dynamics of states with an arbitrary 
number of fermions, including the characteristic interfer- 
ence patterns for indistinguishable fermions in quantum 
mechanics, is correctly described. 

At this point it seems worthwhile to ask some questions 
about the origin of characteristic features of quantum me- 
chanics in our classical statistical setting. Particle-wave 
duality is realized by the discreteness of Ising-spins on one 
side, and the continuous probability distribution or classi- 
cal wave function on the other side. Interference arises from 
the formulation of the basic evolution law in terms of the 
classical wave function. While the classical wave function 
is real, it can nevertheless take positive and negative values 
which can add to zero locally. The superposition principle 
or linearity of the quantum evolution finds a direct origin in 
the formulation of a dynamical law for the classical statis- 
tical ensemble that is linear in the classical wave function. 
The characteristic physics of phases in quantum mechanics 
is connected to the presence of a complex structure within 
the real Grassmann algebra. Planck's constant h appears 
purely as a conversion factor of units. The uncertainty re- 
lations can be obtained directly from the possible solutions 
of the Schrodinger equation. 

Finally, one of the most characteristic elements of the 
mathematical formulation of quantum mechanics is the 
presence of a non-commutative product for operators and 
associated observables. These structures are obviously 
present in our formulation and clearly very useful for a dis- 
cussion of solutions of the Schrodinger equation or Dirac 
equation. We have not addressed in this paper the fun- 
damental origin of non-commuting operators and refer in 
this context to related work [H , [f| , [lj] • An essential point 
is the observation that a classical statistical ensemble ad- 
mits many different product structures for observables, and 
therefore also the definition of many different correlation 
functions. The correct choice of a correlation function for 
the description of a sequence of two measurements depends 
on the details of the measurement process. There exist var- 
ious idealizations of measurements. 

The classical correlation function ( "pointwise multiplica- 
tion" A T B T ) is adapted to a situation of negligible mutual 
influence between two measurements. It is often not appro- 
priate for the description of measurements in subsystems, 
that are characteristic for the quantum experiments. For 
an idealized separation into a subsystem and its environ- 
ment many different classical observables are grouped into 
an equivalence class, for which the properties within the 
subsystem are identical, but the properties relating to the 
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environment differ. Quantum observables are supposed to 
measure only properties of the subsystem. They are there- 
fore associated to a whole equivalence class of classical ob- 
servables, rather than to a single classical observable. A 
one-to-one correspondence between quantum and classical 
observables is a basic assumption of the Kochen-Specker 
theorem Q , which therefore does not apply in our setting. 

The classical product of observables often depends on 
the properties of the environment and is therefore not suit- 
able for idealized measurements of properties of subsys- 
tems. This is the point why Bell's inequalities [H, Q do 
not apply for idealized measurements of properties of sub- 
systems. These inequalities implicitly assume the use of the 
classical correlation function. In short, Bell's inequalities 
are circumvented not by abandoning locality or causality, 
but simply by concentrating on non-classical correlation 
functions that are appropriate for measurements in sub- 
systems. In other words, the experimental verification that 
the observed correlations violate Bell's inequalities tells us 
that the classical correlation function should not be used 
for this type of measurements. This is in accordance with 
our arguments that classical correlations do not correspond 
to idealized measurements of subsystem properties. It has 



been shown that in many circumstances the idealized mea- 
surements of subsystem properties correspond precisely to 
the quantum correlation function that is associated to the 
standard non-commutative operator product in quantum 
mechanics Q, Q. 

The essence of the emergence of non-commutative struc- 
tures is the coarse graining of information. Again, this is- 
sue has not been addressed in the present paper. It seems 
reasonable to expect, however, that a system that is gov- 
erned by the Dirac equation on microphysical scales - say 
lattice distances shorter than the Planck length - will also 
show similar properties at "macroscopic scales" associated 
to coarse graining. (Such macroscopic scales can still be 
much smaller than all characteristic scales of atom physics 
or elementary particle physics.) The basic reason is that 
the form of the Dirac equation for one-particle states is 
essentially fixed by the symmetries. It will not be altered 
if the coarse graining respects the symmetries. While it 
remains an interesting task to perform this coarse graining 
explicitly, the main message of this paper is already very 
clear at the present stage: Quantum field theory can be 
obtained from a classical statistical ensemble. 
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